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ABSTRACT 


Previous  measurements  have  demonstrated  that  a 
polarization  filter  can  increase  ship-background 
temperature  contrast  in  the  infrared,  while  decreasing  the 
received  radiance.  Application  of  this  technique  to 
increasing  range  for  detection  or  recognition  of  ship 
targets  is  being  investigated  through  detection  range 
modeling  for  a  generic  FLIR  sensor.  Laboratory  measurements 
have  been  made  of  effective  Minimum  Resolvable  Temperature 
Difference  (MRTD)  of  a  serial-scan  8-12  |im  sensor  for 
polarized  and  unpolarized  radiation.  A  variety  of  standard 
four-bar  target  boards  of  varied  spatial  frequency  and 
controlled  bar-background  temperature  difference  were  used 
to  construct  MRTD  vs.  spatial  frequency.  Results  were 
compared  with  model  predictions  using  known  or  measured 
component  parameters  for  the  AGA-780  imager,  showing  close 
agreement  for  observations  made  by  a  "trained  observer" .  A 
modified  form  of  MRTD  was  developed  for  a  polarized  target 
using  a  reformulation  of  the  thermal  derivative  of  Planck's 
law.  Modeled  and  measured  values  agreed  closely  for  the 
unpolarized  case,  and  also  for  both  vertically  and 
horizontally  polarized  cases  when  the  appropriate 
parameters  of  the  polarization  filters  were  included. 
Mathematical  analysis  and  measurement  agreed  in  displaying 
an  increase  in  MRTD  with  polarization.  Predictions  of 
maximum  detection  and  recognition  ranges  using  estimates  of 
polarized  effective  target -background  temperature 
difference  indicated  probable  range  improvement  for  sea 
surface  degree  of  polarization  in  excess  of  20%. 
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I .  INTRODUCTION 

With  the  development  of  an  increasingly  sophisticated 
battlefield,  military  systems  have  become  more  dependent  on 
high  technology  than  ever.  Among  other  essential  systems 
necessary  to  achieve  advantage  against  the  opponent, 
thermal  imaging  systems  have  been  demonstrated  to  be 
important  in  numerous  applications.  Detection,  recognition, 
and  identification  from  ground  based,  aerial,  or  space 
platforms  are  among  missions  performed  with  the  use  of 
electro-optic  and  infrared  systems.  Thermal  Imaging  Systems 
rely  upon  a  type  of  signature  found  everywhere  in  nature: 
thermal  radiation . 

Thermal  radiation  exists  in  any  body  at  temperatures 
greater  than  the  absolute  zero  (0°K) ,  which  implies  that, 
not  only  the  target,  but  also  the  background  and  even  the 
atmosphere  will  contribute  to  the  overall  radiation 
reaching  the  sensor.  The  limiting  factor  that  will 
determine  the  possibility  of  performing  some  task  is  known 
as  the  target -to-background  contrast  (TBC) . 

Possible  ways  to  improve  the  target-to-background 
contrast  include  enhancement  of  target  signature  or 
elimination  of  background  radiation.  The  second  one  is  of 
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the  most  interest  for  this  thesis .  Experiments  have 
demonstrated  that,  under  certain  circumstances ,  the  use  of 
polarization  filters  can  reduce  sea  background  [Ref.  1]  . 
The  use  of  polarizers  though,  may  interfere  with  the 
performance  of  the  Thermal  Imaging  System. 

The  most  used  measure  of  performance  for  thermal 
imaging  systems  is  the  Minimum  Resolvable  Temperature 
Difference  (MRTD) ,  which  is  a  laboratory  parameter  defined 
to  determine  the  resolution  of  the  overall  system, 
including  the  operator.  In  order  to  estimate  the  efficiency 
in  performing  a  given  task  under  a  given  scenario,  both  the 
MRTD  and  the  TBC  are  important.  While  the  TBC  is  improved 
when  the  polarizer  is  used,  the  MRTD  values  are  down 
graded.  Trading  off  those  two  parameters  will  determine 
whether  the  probability  of  performing  the  task  is  improved. 

The  objective  of  this  thesis  is  to  investigate  the 
effect  of  the  use  of  polarization  filters  on  the  Minimum 
Resolvable  Temperature  Difference  of  thermal  imaging 
systems.  We  will  start  Chapter  II  by  presenting  some 
fundamentals  of  the  infrared  technology .  Chapter  III  will 
develop  the  mathematical  formulation  of  the  MRTD  for  both 
the  polarized  and  unpolarized  cases .  The  laboratory 
experiment  and  the  analysis  of  the  results  will  be 
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presented  in  Chapters  IV  and  V  respectively.  Chapter  VI 
will  then  summarize  and  conclude  this  work. 
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XI .  FUNDAMENTALS 


The  objective  of  this  Chapter  is  to  ensure  that  the 
reader  understands  the  basic  concepts  of  Infrared 
Radiation.  It  does  not  intend  to  be  a  complete  reference  on 
this  subject,  but  only  to  highlight  those  aspects  of 
Infrared  Radiation  addressed  in  this  thesis. 

A.  ELECTROMAGNETIC  SPECTRUM 

The  electromagnetic  spectrum  is  divided  in  a  "quasi- 
arbitrary"  way.  Its  divisions,  called  bands,  are 
distinguished  from  each  other  mostly  by  the  methods  used  to 
produce  and  detect  the  radiation. 

The  infrared  portion  of  the  spectrum  lies  in  the 
region  between  the  visible  and  the  microwave,  and  is 
subdivided  in  four  regions.  The  first,  from  0.7|im  to  1.1pm, 
is  the  near  infrared  (NIR)  and  is  dominated  by  the 
reflected  sun  radiation.  Low  light  level  TVs,  image 
intensifiers,  and  night  vision  devices  operate  in  this 
region  [Ref.  2].  The  second  region  spans  from  1 .  l|im  to  2.5|Hm 
and  is  called  the  short  wavelength  infrared  band  (SWIR) . 
The  third  region,  the  mid-wavelength  infrared  (MWIR) , 
ranges  from  2.5pm  to  7.0pm.  It  is  common  to  consider  this 
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region  as  spanning  from  3  pm  to  5pm,  since  atmospheric 
attenuation  limits  its  useful  range.  The  next  region  (LWIR) 
is  the  most  important  for  this  thesis,  since  it  is  where 
all  the  laboratory  measurements  were  taken.  It  ranges  from 
7pm  to  15pm,  but  is  also  limited  from  8pm  to  12pm  for 
practical  use. 

Figure  2.1  shows  a  pictorial  representation  of  the 
electromagnetic  spectrum,  with  special  emphasis  on  the  IR 
spectrum. 
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Figure  2.1  -  The  Electromagnetic  Spectrum  (After  Ref.  5) 


It  is  important  to  emphasize  that,  being  part  of  the 
electromagnetic  spectrum  means  that  the  infrared  radiation 
obeys  all  the  laws  of  electromagnetism,  as  stated  by  the 
Maxwell's  Equations: 
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(2.1) 


AxE=~ir 

dt 


AxH  =  J  + 


A  •  D  =  pv 


dD_ 

dt 


(2.2) 

(2.3) 


A*B  =  0 


(2.4) 


where : 

E  =  electric  field  intensity  vector  (V/m) 
H  =  magnetic  field  intensity  vector  (A/m) 
D  =  electric  flux  density  vector  (C/m2) 

B  =  magnetic  flux  density  vector  (T) 

J  =  electric  current  density  (A/m2) 
t  =  time  (s) 

pv  =  volume  charge  density  (C/m3) 


B.  LAWS  OF  THERMAL  RADIATION 

Before  starting  to  enumerate  the  main  laws  of  the 
thermal  radiation,  some  definitions  should  be  presented: 

•  Emissivity  (8)  :  it  is  "the  ratio  of  the  radiant 
exitance  or  radiance  of  a  given  body  to  that  of  a 
blackbody."  [Ref.  3] 
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•  Blackbody :  "it  is  an  object  which  absorbs  all 
radiation  that  impinges  upon  it  at  any  wavelength." 
Blackbodies  have,  by  definition,  emissivity  equal  to 
unity  (£  =  1)  [Ref.  4]. 

•  Gray  body:  an  object  whose  emissivity  is  high  and 
fairly  constant  with  wavelength.  [Ref.  5] 

Table  2.1  shows  the  most  commonly  used  infrared 
quantities . 


Quantity 

Symbol 

Units 

Definition 

Radiant  Energy 

U 

Joules 

Radiant  Power 

P 

Watts 

8u/5t 

Radiant  Emittance 

w 

Watt /cm2 

8p/3s 

Radiant  Intensity 

J 

Watt/sr 

8p/3Q 

Radiance 

N 

Watt /cm2  sr 

d2P/  (dSdQ.) 

Table  2.1  -  Radiometric  Quantities  (After  Ref.  6). 


1.  Planck's  Law 

It  describes  the  spectral  distribution  of  the 
radiation  from  a  blackbody,  and  is  given  by  the  following 
formula : 


lithe 

X\ehc'm-\) 


xlO 


,-6 


(2.5) 
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where : 


WjL,b  =  the  blackbody  spectral  radiant  emittance  at 
wavelength  X  (Watt/m2  jim) 

c  =  the  velocity  of  light  (3  x  108  m/s) 
h  =  Planck's  constant  (6.6  x  10"34  Joule/s) 
k  =  Boltzmann's  constant  (1.4  x  10~23  Joule/K) 

T  =  absolute  temperature  of  the  blackbody  (K) 

X  =  wavelength  (m) 

2.  Wien's  Displacement  Law 

By  differentiating  Planck's  formula  with  respect  to 
the  wavelength  (X)  ,  and  finding  the  maximum,  we  can 
establish  the  wavelength  where  the  peak  of  radiation  occurs 
for  a  given  temperature.  Wien's  law  is  stated  as: 


2898 

T 


(2.6) 


where : 

Xroax  =  wavelength  where  the  peak  of  radiation  occurs 

(Mm)  . 

This  formulation  can  be  graphically  observed  in  Figure 
2.2,  where  both  Planck's  and  Wien's  laws  are  plotted. 
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Wavelength  -  micrometer 


Figure  2.2  -  Spectral  Radiant  Exitance  of  Blackbodies  at 

Various  Temperatures. 


3 .  Stef an-Boltzmann  Law 

We  can  obtain  the  total  radiant  emittance  of  a 
blackbody  by  integrating  Planck's  law  over  the  whole 
spectrum  (from  X=0  to  X=°°)  .  It  is  known  as  the  Stefan- 
Boltzmann  Law,  and  is  formulated  as: 

Wb=(fT 4  (2.7) 

where : 

Wb  =  total  radiant  emittance  of  a  blackbody 

<7  =  the  Stef  an-Boltzmann  constant  (5.7  x  10”8  Watt/m2) 
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4. 


Total  Power  Law 


"When  radiation  is  incident  upon  a  body,  some  of  it  is 
transmitted,  some  absorbed,  and  some  is  reflected.  Thus  the 
ratios  of  each  of  these  to  the  incident  power  must  add  up 
to  unity."  [Ref.  3] 

a+p+z  =  \  (2.8) 

where : 

a  =  absorptivity  =  Pabsorbed/Pincident 
P  =  reflectivity  =  Preflected/Pincident 

T  =  transmissivity  =  Ptransmitted/Pincident 

5.  Kirchof f ' s  Law 

"Kirchoff's  law  of  electromagnetic  radiation  states 
that  a  good  absorber  is  also  a  good  emitter  of  radiation  by 
the  equation:"  [Ref.  6] 


C.  ATMOSPHERIC  PROPAGATION 

Four  major  phenomena  occur  when  electromagnetic 
radiation  propagates  from  a  target  to  a  sensor  through  the 
atmosphere:  (1)  its  intensity  is  reduced,  (2)  atmospheric 
radiance  is  added  to  the  target  radiance,  (3)  non-scene 
radiance  is  scattered  into  the  sensor,  and  (4)  some  scene 
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radiance  is  scattered  away  from  the  sensor .  Those  phenomena 


are  illustrated  in  Figure  2.3. 


Field  of  View 

Figure  2.3  -  Atmospheric  Effects  (After  Ref.  2). 

Two  values  are  commonly  used  to  define  the 
atmosphere's  behavior.  The  extinction  coefficient  (y(A.)  ) 
represents  the  total  reduction  of  radiation  in  the  path 
from  the  source  to  the  sensor,  whereas  transmittance  (x(A,)) 
is  defined  as  the  ratio  between  the  received  radiance  at 
the  sensor  to  the  emitted  radiance  from  the  source.  It  is 
not  difficult  to  measure  those  values  but  it  is  very 
cumbersome  to  predict  them,  since  the  atmosphere  is  very 
dynamic  and  complex.  Models  have  been  proposed  and  some 
computer  codes  such  as  LOWTRAN,  MODTRAN,  HITRAN1,  and 
SEARAD2  are  used  to  obtain  y(X)  or  t(?i)  . 


1  LOWTRAN,  MODTRAN,  and  HITRAN  -  Low,  Moderate,  and  High  Resolution  Transmission  Code 
developed  by  the  US  Air  Force  (Geophysics  Laboratory)  . 

2  SEARAD  -  Modification  of  MODTRAN  developed  by  the  US  Navy  (NRaD)  . 
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In  this  thesis,  SEARAD  will  be  used  to  predict 
atmospheric  transmittance  whenever  needed.  Figure  2.4  shows 
an  example  of  transmittance  values  obtained  from  SEARAD.  We 
can  observe  that  wavelength  severely  affects  the 
transmittance  through  the  atmosphere,  leaving  us  with  two 
practical  windows  to  work  in:  3-5fim  (MWIR)  and  8-12|Lim 
(LWIR)  . 


wavelength  -  micrometer 


Figure  2.4  -  Atmospheric  Propagation  (1  Km  Horizontal  Path 
Length,  Midlatitude -Summer,  Navy  Maritime  Aerosol,  Airmass 
Characteristic  =3,  Calculated  Using  Searad) . 


D.  POLARIZATION 

Electromagnetic  waves  (such  as  Infrared  radiation)  are 
composed  of  an  electric  and  a  magnetic  field.  Those  fields 
are  perpendicular  to  each  other  and  perpendicular  to  the 
direction  of  propagation.  Polarization  of  the  wave  is 
defined  as  the  position  of  the  electric  vector  as  it 
progresses  [Ref.  7]. 
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The  polarization  of  the  wave  may  change  in  the 
boundary  between  two  media.  This  phenomenon  was  studied  by 
Fresnel,  who  proposed  the  following  equation  for 
nonmagnetic  media: 


sin  (0,-0,) 

V  .  ,  r= - — 

(2  . 

.10) 

perpendicular  ^(0.+^) 

tan(0(  —0t) 

(2. 

.11) 

Wi"  tan(0.+0;) 

2  •  sin(0, )  •  cos(0; ) 

(2. 

.12) 

^perpendicular  Sjn(^.+^) 

2-sin(0t)-cos(0,) 

(2. 

.13) 

*****  ~  sin (6i  +6t)  cos(0i  -dt) 

where : 

0i  =  incident  angle. 

0t  =  angle  of  transmission, 
r  =  reflection  coefficient, 
t  =  transmission  coefficient. 

We  must  also  define  degree  of  polarization  as: 


POL  = 


<*,)+{*„} 


(2.14) 


where : 

<NV>  =  mean  apparent  vertically  polarized  radiance. 
<Nh>  =  mean  apparent  horizontally  polarized  radiance. 
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The  Fresnel  Equations  demonstrate  that,  after  reaching 
the  boundary  between  two  media,  unpolarized  incident 
radiation  may  be  reflected  with  a  certain  degree  of 
polarization,  since  the  amount  of  parallel  and 
perpendicular  polarization  components  are  reflected  obeying 
different  equations.  At  some  angle  (called  the  Brewster 
angle)  ,  we  may  have  rparanei  equal  to  zero  and  rperpendicuiar 
greater  than  zero.  In  this  case,  unpolarized  radiation 
would  reflect  vertically  polarized,  which  means  a  degree  of 
polarization  greater  than  zero. 
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III.  MINIMUM  RESOLVABLE  TEMPERATURE  DIFFERENCE 


Minimum  resolvable  temperature  difference  (MRTD)  is 
defined  as  the  temperature  difference  between  the 
background  and  a  set  of  four  standard  bars  (7:1  aspect 
ratio)  required  to  make  the  bars  just  resolvable  [Ref.  8]. 
It  is  the  most  used  measure  of  performance  for  thermal 
imaging  systems . 

Before  giving  its  complete  mathematical  formulation, 
it  is  necessary  to  highlight  some  important  characteristics 
of  the  MRTD.  First,  it  is  a  function  of  the  entire  system, 
including  the  observer.  This  implies  the  inclusion  of  some 
psychophysical  aspects  of  human  vision  in  the  final 
formulation  of  the  MRTD.  Another  important  feature  of  the 
MRTD  is  that  it  is  a  function  of  the  spatial  frequency 
(Figure  3.1).  It  allows  us  to  predict  field  performance  of 
the  thermal  imaging  system  by  establishing  some  criteria  to 
convert  object  dimensions  into  spatial  frequency. 

There  are  several  formulations  used  for  the  MRTD.  In 
this  work,  we  will  use  an  approach  similar  to  that  proposed 
by  Lloyd  [Ref.  9]  .  Some  concepts  such  as  the  signal-to- 
noise  ratio  (SNR)  and  the  noise  equivalent  temperature 
difference  (NETD)  will  be  introduced  prior  to  the  MRTD. 
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Spatial  fT  2fT  3fT  4fT  5fi 

frequency 

Figure  3.1  -  MRTD  Patterns  of  Differing  Spatial  Frequency. 


A.  SIGNAL-TO-NOISE  RATIO  AND  NOISE  EQUIVALENT  TEMPERATURE 


DIFFERENCE 


The  first  concept  to  be  developed  is  the  importance  of 
the  signal-to-noise  ratio  (SNR)  in  the  thermal  imaging 
process.  The  SNR  gives  a  fundamental  limit  to  the 
performance  of  a  thermal  imaging  system.  SNR  is  defined  as 
the  relationship  between  the  signal  and  the  noise  level  at 
the  output  of  the  system.  It  interferes  with  the  observer's 
perception  of  the  scene.  It  is  usually  used  in  an  inverse 
form,  as  the  noise  equivalent  temperature  difference 
(NETD) .  NETD  is  the  temperature  difference  between  the 
signal  and  the  background,  which  produces  a  SNR  of  one  in 
the  video  signal. 


NETD  = 


AT  _  AT 
VjVn  "  SNR 


(3.1] 
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where : 


NETD  =  Noise  Equivalent  Temperature  Difference  (°K) 

AT=  Target-Background  Temperature  Difference  (°K) 

Vs  =  rms  output  signal  voltage  (V) 

Vn  =  rms  output  noise  voltage  (V) 

SNR  =  Signal-to-Noise  Ratio 

Noise  equivalent  temperature  difference  is  obtained  by 
considering  a  large  target,  surrounded  by  a  large 
background  (Figure  3.2),  both  with  uniform  temperature  (Tt 
and  Tt,  respectively)  . 

Tb 

I  Tt  I 


radiance 

(3.2) 


Figure  3.2  Target  and  Background  Geometry. 

Assuming  a  Lambertian  surface,  the  spectral 
of  the  surface  is: 


where : 

Njl  =  spectral  radiance  (watt/cm2  |im  sr) 

=  spectral  radiant  exitance  (watt/cm2  |im) 
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The  instantaneous  field  of  view  subtended  by  the 
detector  area  (axb)  is  axp,  so  that  the  detector  receives 
power  from  the  area  otxpxR2  in  .  the  object  plane.  The 
spectral  power  at  the  aperture  (Figure  3.3)  is  then  given 
by  the  expression: 


W. 


WA?P£l 


aperture 


nR 


w. 


WAficfi 


aperture 


K 


(3.3a) 


(3.3b) 


where : 

Waperture  =  spectral  power  at  the  aperture 

A0  =  aperture  area  (m2) 

a,  b  =  detector  dimensions  (m) 

a,  P  =  detector  instantaneous  field  of  view  (rad) 


20 


Therefore,  the  power  received  by  the  detector  is: 


f,=— A'M-  *•„(•*) 

7t 


(3.4) 


where : 

Px.  =  power  received  by  the  detector  (watt/|im) 
x0  =  transmissivity  of  the  optical  system 


Hence,  the  differential  power  is: 


^  =  M  4  ,  M 

ar  n  0  0  dT 


(3.5) 


which  gives  a  signal  voltage  of: 


ar  k 


K-*o 


dT 


■R 


(3.6) 
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where  R  (Responsivity  in  V/W)  is  given  by: 


v.-p'q) 

A/T 


(3.7) 


Vn  =  noise  voltage  in  the  bandwidth  Afn(V) 
D*  =  Detectivity  (cm  Hz%  W"1) 


Integrating  Equation  3.4  over  the  wavelength,  we  have: 


V,  W.  i  sr 


(3.8) 


Since  we  defined  NETD  as  the  AT  for  a  SNR  of  one. 


then: 


NETD  =  - 


nJa  ■  b  ■  Afn 


cc-p-A0-]^D'a)-T0{X)-ca 


(3.9) 


For  a  system  filtered  to  pass  a  selected  band  (Aa  to 
X.2 )  ,  and  assuming  xQ  constant  within  that  band. 


Uwa 


-r0(X) -dX  =  tb  dX 


(3.10) 
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We  define  effective 


dW 

dT 


so  that 


AW 

AT 


4 


■D\Xp)  =  D\Xp)-j 


dW  D*(X) 
dT  D*(Xp) 


dX 


(3.11) 


where : 

D*(Ap)  =  detectivity  at  the  peak  wavelength  (cm  Hz'*  W"1) 
and : 


NETD  = 


a-P-A0-t0-D\Xp) 


AW 

AT 


(3.12) 


Now,  the  Universal  Radiation  Curve  can  be  used  to 


evaluate 


AW 

AT 


,  and  also, 


CO 

6 


(3.13) 


In  Equation  3.13,  Afn  is  the  noise  equivalent  bandwidth 
of  the  reference  filter,  or  equivalent  external  filter 
including  the  effect  of  all  electronic  MTF's,  S(f)  is  the 
noise  spectrum  (normalized)  ,  and  H(f)  is  the  Modulation 
Transfer  Function  of  the  reference  filter.  This  is  related 
to  the  sensor  dwell  time  ta- 
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#(/)= 


(3.14) 


with 


f  MB  ~ 


2-tj 


(3.15) 


The  filter  is  then  suitably  limited  for  the  scanning 
rate,  where  twice  the  dwell  time  is  equivalent  to  one 
period.  Assuming  the  system  has  n  detectors, 


(3.16) 

FOVx  ■  FOVy  •  Fr 


where : 

n  =  number  of  detectors 
T|sc  =  scan  efficiency 

FOVx  and  FOVy  =  overall  field  of  view  (mrad) 

Fr  =  frame  rate  (Hz) 

We  can  finally  write  a  formula  for  NETD  that  is 
convenient  for  trade-off  analysis: 


NETD 


20  ■  ^7T  a  b-  FOVx  •  FOVy  ■  Fr 

AW 


AT 


fvM-f.’o’a,)''. 


(3.17) 
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where : 


D0  =  aperture  diameter  (m)  is  equal  to 

B.  MINIMUM  RESOLVABLE  TEMPERATURE  DIFFERENCE  (MRTD) 

As  defined  at  the  beginning  of  this  Chapter,  MRTD  is  a 
measure  of  the  observer's  ability  to  discern,  at  the  output 
of  a  thermal  imaging  system  (usually  a  gray  scale  display) , 
a  given  spatial  frequency.  It  includes  some  psychophysics 
of  the  human  eye,  such  as  eye/brain  temporal  integration 
and  the  effective  narrow  band-spatial  filter  of  the  eye. 

The  following  assumptions  are  made: 

•  Spatial  filtering  in  the  eye,  for  a  given  target 
spatial  frequency  (fT)  ,  is  approximated  by  a  matched 
filter  with  the  following  transfer  function  [Ref.  9, 
p.183]  : 

H  f  =  sine  —  (3.18) 

2  fT 

•  The  amplitude  of  the  fundamental  frequency  in  the 
square  wave  is  A/n  times  the  square  wave  amplitude, 

*„=I[a«7vJ  (3-19) 

K 
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where : 


Rsw  =  square  wave  response  of  the  system 

MTFsystem  =  modulation  transfer  function  of  the  whole 

system 

When  the  observer  is  introduced,  some  perception 
factors  should  be  included  [Ref.  9]: 

•  Eye  signal  averaging:  the  eye  senses  the  mean  signal 
of  the  image,  which  is  2/ K  for  a  half -cycle  of  the 
sine  wave. '  This  factor  associated  with  the 
assumption  stated  in  Equation  3.19  leads  to  a  factor 
of  8/7C2  when  estimating  the  square  wave  by  using  a 
sine  wave; 

•  eye  integration  time  -  Eye/brain  integration  time 
(Te)  is  assumed  to  be  approximately  0.2s,  which 
improves  the  SNR  by  a  factor  of  (TeFr)I/4  [Ref.  9]  . 

•  the  eye  integrates  the  image  along  the  bar  (Figure 
3.4).  It  gives  an  improvement  factor  in  the  signal- 
to-noise  ratio  of: 

(lTJVwT  f  7  r  (3.20) 

UJ  "U  J  i2-/r-*J 
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Figure  3.4  -  Eye  Integration  along  the  Bar. 


•  eye  spatial  filtering:  in  order  to  account  for 
the  monitor  and  the  observer,  the  new  noise 
equivalent  bandwidth  becomes 


A/„'  =  JS(/)  •  •  MTF^, 

0 


f 

•sine 

v 


■df 


(3.21) 


which  corresponds  to  a  SNR  improvement  of : 


P 


1/2 


(3.22) 
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Now,  if  we  let  the  SNRT  be  a  minimum  signal-to-noise 
ratio  needed  for  detection  of  an  individual  bar,  we  can 
solve  for  AT  as 


NETD-p 


1/2 


a  r=- 


(hM} 

7 


1/2 


■snrt 


(3.23) 


t  jO 

y7T2  j 


MTF^fcFrJ 


/  2 


The  probability  of  detecting  a  single  bar  is  a 
function  of  the  SNR.  Then,  calculating  the  temperature 
difference  corresponding  to  a  desired  SNR  gives  us  the 
Minimum  Resolvable  Temperature  Difference  (MRTD) 


MRTD(f)  = 


NETD-p112  if Tfi)112  - SNR 
1.52  ■MTFmm  (T,Fr)'n  ' 


(3.24) 


Finally,  for  white  noise,  pA  can  be  approximated  as 


,1/2 


2 Jr 

V4 fn 


(3.25) 


and  the  final  expression  for  the  MRTD  becomes 


MRTD(fr)  = 


NETD  •  SNRT  ■  fT  {cc0fn  • 

1 .52  •  MTFsystem  ■  (TeFrTdAfl if'2' 


(3.26) 
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which  is  equivalent  to  Equation  5.58  in  Lloyd  [Ref.  9]. 

Various  authors  have  proposed  some  alternative 
expressions  for  the  MRTD.  Shumaker  et  al  [Ref.  8]  give  the 
following  formula: 


MRTD(f)  = 


20  SNRt  •  (FOVx  •  FOVy)U2  f  •  pf_ 


to  D-D* -{it -n  TjJ2  •  ■  MTF^  ■{. LTe  )U2 


3N 


(3.27) 


where : 

SNRt  =  perceived  signal-to-noise  ratio  set  according  to 

the  desired  probability  of  performing  the  given  task 

f  =  spatial  frequency  (cycles/mrad) 

D  =  aperture  diameter  (m) 

D*  =  detectivity  (cm  Hz%  W_1) 

n  =  number  of  detectors 

T|sc  =  scan  efficiency 

x0  =  transmission  of  the  optics 

a  =  in-scan  detector  angular  subtense  (mrad) 

P  =  cross-scan  detector  angular  subtense  (mrad) 

MTF system  =  the  system  modulation  transfer  function 

L  =  the  length- to-width  ratio  for  the  bar  chart  (7) 

Te  =  eye  integration  time  (0.2s) 
dN 

=  derivative  of  the  Planck's  Law  (watt  cm-2  K'1  sr"1) 
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FOVx  =  in-scan  field  of  view  (mrad) 

FOVy  =  cross-scan  field  of  view  (mrad) 
px  =  noise  filter  factor 

Equation  3.27  corresponds  to  substituting  the 
expression  for  the  NETD  (Equation  3.17)  in  the  formula  for 
the  MRTD  (Equation  3.26). 

In  the  next  section,  we  will  study  the  effects  of  the 
inclusion  of  a  polarizer  in  the  system. 

C.  MRTD  WITH  THE  POLARIZER 

In  this  section,  we  will  propose  a  formulation  for  the 
MRTD  for  the  case  when  a  polarizer  is  added  to  the  thermal 
imaging  system.  For  simplicity,  from  now  on,  we  will  call 
it  MRTDp .  In  this  study,  all  the  parameters  involved  in  the 
MRTD  calculation  were  individually  analyzed  and  their 
potential  to  affect  the  value  of  the  MRTDP  evaluated.  The 
parameters  considered  to  represent  the  greatest  impact  in 
this  formulation  are  the  polarizer  transmissivity  (Tp) , 
polarizer  Modulation  Transfer  Function  (MTFP)  ,  and  the 
partial  derivative  of  Planck's  Equation  in  relation  to  the 
temperature.  They  will  be  discussed  in  the  rest  of  this 
section. 
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1.  "Thermal  Derivative" 

The  "thermal  derivative"  is  the  partial  derivative  of 
Planck's  law  with  respect  to  temperature.  It  is  a  function 
of  wavelength  and  the  background  temperature.  Thermal 
imaging  systems,  in  fact,  do  not  measure  temperature  but 
radiance  differences.  Temperature  difference  (AT)  is, 
however,  a  convenient  concept  and  can  be  used  when  both  the 
system  spectral  response  and  the  background  temperature  are 
specified.  [Ref.  2] 

The  addition  of  a  polarizer  in  the  system  implies  that 
a  fraction  of  the  total  radiance  will  not  reach  the 
detector.  For  a  perfectly  unpolarized  image  (target  and 
background) ,  and  assuming  a  perfect  polarizer  (efficiency 
equal  to  unity),  this  reduction  would  be  fifty  percent. 
However,  in  some  situations,  the  concept  of  a  perfectly 
unpolarized  image  is  not  met.  For  example,  "Strong  vertical 
polarization  is  observed  in  the  sea  surface  emission  near 
the  Brewster  angle  in  the  8-12|jm  (LWIR)  band,  whereas 
horizontal  polarization  due  to  reflection  is  usually 
dominant  in  the  3-5pm  (MWIR)  band."  [Ref.  1]  Exploring  this 
effect  would  allow  us  to  improve  the  overall  contrast  by 
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using  a  polarizer  that  filtered  out  the  unwanted  background 
radiation. 

Since  the  "thermal  derivative"  is  usually  calculated 
and  tabulated  for  a  standard  background  temperature  of  300 
degrees  Kelvin,  our  task  is  to  estimate  the  value  of 

background  temperature  that  would  reduce  the  total  in-band 
radiance  to  one  half  of  its  previous  value.  In  other  words, 
we  would  need  to  solve  the  following  equation 

rdN,(T B\dX=k-  [ (3.28) 
{  {  3T 

where : 

Tb  =  background  temperature  (K) 

k  =  factor  that  represents  the  fraction  of  the  total 

in-band  radiance  that  reaches  the  detector,  when  the 

polarizer  is  used  (0.5  for  a  perfect  polarizer  and  a 
perfectly  unpolarized  image).  When  k  is  equal  to  0.5,  the 
equivalent  value  of  TB  is  283K  for  the  3 -5pm  band  and  262K 

for  the  8-12|im  band.  (Ref.  8] 

2 .  Polarizer  Transmissivity 

Another  impact  of  the  polarizer  on  the  overall 
calculation  is  due  to  the  fact  that  its  transmissivity  is 
different  from  unity.  In  fact,  the  transmissivity  of  the 
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polarizer  used  in  this  experiment  (aluminum  grid 
superimposed  on  KRS  substrate  manufactured  by  Graseby- 
Specac)  is  given  by  the  manufacturer  as  being  0.85,  for  the 
electric  vector  aligned  with  the  passage  direction. 

Considering  that  the  impact  of  the  loss  in 
transmissivity  would  be  equivalent  to  degrading  the  overall 
transmissivity  of  the  optics,  we  would  have 


MRTDp 


MRTD 


(3.29) 


where : 


Xp  =  polarizer  transmissivity. 

3 .  Polarizer  MTF 

Assuming  the  thermal  imaging  system  is  linear,  its 
modulation  transfer  function  is  the  result  of  the 
multiplication  of  the  MTF  of  the  individual  subsystems 
(optics,  stabilization  element,  detector,  amplifiers, 
electronic  filters,  LEDs,  visual  optics  and  eye)  [Ref.  2]. 
The  addition  of  the  polarizer  to  the  system  will  include 
one  more  factor  to  the  MTF 


MTF 


polarized  system 


=  MTFr 


polarizer 


■MTF 


unpolarized  system 


(3.30) 


and  it  would  affect  the  MRTD. 
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MRTD 


MRTDp 


MTFpolarizer-tp-k 


(3.31) 


Since  the  equipment  necessary  to  measure  the  MTF  of 
the  polarizer  was  not  available,  we  estimated  it  by 
measuring  the  MTF  of  the  overall  system  with  and  without 
the  polarizer.  The  results  are  shown  in  Chapter  IV. 

Since  all  values  in  the  denominator  of  Equation  3.31 
are  less  than  one,  we  can  expect  the  curve  of  the  polarized 
MRTD  to  be  shifted  upward  in  relation  to  the  unpolarized 
MRTD.  This  effect  will  be  observed  in  the  next  chapter  when 
we  calculate  the  MRTD  and  the  MRTDP  for  the  AGA  780. 


34 


IV 


LABORATORY  EXPERIMENTS 


A.  AGA  780 

The  AGA  Thermovision®  780  is  a  fourth  generation  dual¬ 
scanner  thermal  imaging  system  manufactured  by  AGA  Infrared 
Systems  AB,  with  headquarters  in  Danderyd,  Sweden.  In  this 
experiment,  only  the  LWIR  band  (8-12|im)  was  used  and  the 
system  characteristics,  when  available,  were  obtained  in 
the  operating  manual.  Some  parameters  were  not  available 
and  had  to  be  estimated  using  typical  values  obtained  in 
the  literature. 

B.  CALIBRATION 

As  mentioned  in  Chapter  I,  the  rate  of  emission  of 
infrared  radiation  from  a  blackbody  increases  non-linearly 
with  temperature.  It  is  necessary  to  establish  the  correct 
relationship  between  actual  object  temperature  and  thermal 
imaging  output.  This  is  done  by  using  calibration  curves. 

In  the  operation  of  the  AGA  780,  objects  are 
characterized  by  temperature  differences,  but  the  analog 
output  of  the  system  is  given  in  "isothermal  units." 
Calibration  curves  are  then  used  to  transform  isotherm 
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units  into  temperature  for  each  configuration  of  the 
equipment  (aperture  and  lenses) .  When  digitization  boards 
and  digital  data  acquisition  software  are  used,  the 
fundamental  output  units  will  be  "digital  levels . " 

Digital  Level  (DL)  is  the  scale  used  by  the  signal 
processing  computer,  and  is  dependent  on  the  ADC  boards  and 
software  used.  For  the  PTRWin  [Ref.  17]  Data  Acquisition 
and  Processing  System  currently  used  with  the  AGA  780 
imager  at  the  Naval  Postgraduate  School  this  is  based  on  a 
12-bit  output,  ranging  from  0  to  4095  (212  =  4096  levels)  . 

Digital  level  values  can  be  converted  to  Isothermal  Units 
(IU) .  Isothermal  units  depend  on  two  settings  of  the 
thermal  imaging  system,  the  thermal  range,  and  the  thermal 
level,  and  range  from  zero  to  2000.  By  setting  a  thermal 
level,  we  establish  the  midscale  level  of  thermal  units 
allowed  in  the  output.  On  the  other  hand,  by  setting  the 
thermal  range,  we  establish  the  variation  from  the  midscale 
allowed  in  the  output.  For  example,  setting  the  thermal 
level  to  100  and  the  thermal  range  to  20,  we  will  have  the 
output  ranging  from  90  to  110  IU.  The  thermal  range  is  then 
divided  into  discrete  steps,  as  allowed  by  the  available 
digital  levels.  In  the  previous  example,  20  IU  is  divided 
by  4095,  to  establish  the  correspondence  that  1DL  is 
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equivalent  to  0.00488IU.  The  relationship  between  Digital 
Level  and  Isothermal  Units  is  given  by: 


TTT  DL- Thermal  Range  _  ,  _  .  Thermal  Range 

IU  = - = - —  +  Thermal  Level - - - — 

4095  “  2 


(4.1) 


where : 

IU  =  Isothermal  Units 

DL  =  Digital  Level. 

The  conversion  between  Isothermal  Units  and 
Temperature  in  now  obtained  through  calibration  curves.  In 
order  to  establish  the  calibration  curves,  laboratory 
measurements  were  performed  following  the  steps  given  in 
the  AGA  780  Operating  Manual.  Basically,  a  number  of 
blackbody  sources  at  different  temperatures  were  measured 
with  the  system.  The  results  were  fed  into  a  Mat  lab® 
program  (Appendix  A)  that  gives  the  least  square  fitting 
parameters  for  the  following  expression  [Ref.  4] 


IU  = 


A 

B 


C-eT  -1 


where : 

IU  =  Isothermal  Units 
T  =  Absolute  Temperature 


(4.2) 
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A,  B,  C 


Calibration  constants  that  depend  on  actual 


aperture,  filter  scanner  version,  etc. 

Table  4.1  shows  the  results  obtained  from  these 
measurements,  and  the  corresponding  values  of  A,  B,  and  C. 

A  Matlab®  function  was  written  to  perform  conversions 
between  Digital  Level  and  Temperature.  This  function, 
dl2temp.m,  is  shown  in  Appendix  A. 


Temperature 

(°C) 

Filter 

Temperature 

(°C) 

Filter 

Unpolarized 

Vertical 

Polarized 

Horizontal 

Polarized 

22 

51 

23 

50 

50 

27 

55 

27 

51 

51 

32 

59 

33 

52 

52 

37 

63 

37 

54 

54 

42 

68 

43 

57 

57 

50 

75 

50 

59 

59 

60 

84 

61 

63 

63 

70 

94 

70 

67 

67 

80 

103 

80 

71 

71 

91 

115 

90 

75 

75 

101 

126 

100 

80 

80 

111 

139 

110 

86 

86 

A 

2 . 3 18E+3 

A 

5.679  E+3 

5.679  E+3 

B 

1.218E+3 

B 

57.46 

57.46 

C 

0.742 

C 

0.918 

0.918 

Table  4.1  -  Output  in  Isothermal  Units. 


C.  MODULATION  TRANSFER  FUNCTION 


As  observed  in  Equations  3.27  and  3.32,  the  Modulation 
Transfer  Function  (MTF)  greatly  impacts  the  final  value  of 
the  MRTD .  One  important  aspect  to  be  pointed  out  is  the 
fact  that  the  MTF  goes  to  zero  at  some  value  called  the 
cutoff  spatial  frequency  (Figure  4.4).  It  implies  that  the 
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MRTD  curve  will  tend  to  asymptote  at  this  cutoff  frequency 
(Figure  4.5).  In  practice,  most  current  thermal  imaging 
systems  have  this  cutoff  at  a  frequency  near  the  reciprocal 
Of  the  detector  angular  subtense  [Ref.  8] .  The  next  section 
will  show  how  the  cutoff  frequency  for  the  AGA  780  was 

experimentally  estimated. 

1.  Cutoff 

During  the  preliminary  measurements,  it  was  observed 
that  spatial  frequencies  greater  than  0 . 6  cycles  per 

milliradian  could  not  be  resolved  by  the  system.  In  order 
to  determine  the  exact  value  of  the  cutoff,  the  following 
approach  was  adopted.  At  2  meters  and  using  a  standard 

four-bar  pattern  that  had  a  bar  width  of  5mm,  the  output 

was  recorded.  Figure  4.1  shows  the  recorded  profile  of  this 
output,  using  the  "profile  function"  of  the  PTRWin® 
software  [Ref.  17] . 


Figure  4.1  -  Output  of  the  System  at  0.4cycles/mrad. 

This  corresponds  to  an  angular  resolution  of 

5mm/2m  =  2.5  mrad,  which  would  give  a  cutoff  of 
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0.4  cycles /mrad.  Since  the  first  peak  still  has  a  short 
flat  top,  we  would  infer  that,  in  fact,  the  minimum  angular 
resolution  would  be  greater  then  0.4  cycles/mrad.  Then,  at 
2.75m  and  using  a  pattern  with  a  spatial  frequency  of 
0.55  cycles/mrad),  we  had  the  profile  shown  in  Figure  4.2. 
This  figure  shows  that  a  5mm  bar  is  represented  by  a  single 
point,  which  indicates  that  the  cutoff  frequency  of  the 
system  is,  indeed,  approximately  0.55  cycles/mrad. 


Figure  4.2  -  Output  of  the  System  at  0.55  cycles/mrad. 
2.  MTF 


The  next  step  necessary  to  calculate  the  MRTD  of  the 
system  was  to  determine  the  Modulation  Transfer  Function 
(MTF)  .  MTF  is  a  function  of  the  spatial  frequency  that 
defines  the  ratio  of  the  modulation  in  the  observed  image 


to  that 

in  the  actual 

object,  and  is 

defined 

as 

in 

Equations 

4 . 3  and  4 . 4 

[Ref.  5] .  Figure 

4 . 3  shows 

the 

difference 

in  modulation 

between  the  input 

and  the 

output 

the  system. 


M(/)  = 


V  -V  ■ 

T  max  min 

V  +V  ■ 

r  max  r  min 


(4.3) 
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MTF  -  Mpuif) 

MJf) 


(4.4) 


where : 

M(f)  =  Modulation 

Min(f)  =  Modulation  at  the  input 

Mout(f)  =  Modulation  at  the  output 

f  =  spatial  frequency 

Vmax  and  Vmin  as  shown  in  Figure  4.3. 


-10  -5  0  5  10 

units  of  length 

(a) 


-10  -5  0  5  10 

units  of  length 

(b) 


Figure  4.3  -  Output  Signal (b)  Resultant  of  the  Impact  of  the 

MTF  on  the  Input  Signal (a). 

For  the  purpose  of  this  thesis,  the  MTF  of  the  AGA  780 
was  estimated  by  measuring  and  recording  the  output  of  the 
system  at  several  selected  spatial  frequencies  and 
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calculating  the  modulation  of  each  one .  The  result  is  shown 
in  Figure  4.4.  The  following  step  was  to  find  a  curve  that 
best  fitted  the  data.  Using  Matlab®,  we  generated  a  fitting 
algorithm  (mtffit.m.  Appendix  A)  that  tried  the  several 
templates  for  the  fitting  curve: 

•  polynomial  of  degree  2  or  3; 

I 

•  gaussian; 

•  typical  lowpass  filter;  and 

•  "sine"  function. 

Using  the  least  square  approach,  the  template  that 
best  fitted  the  data  was 


( .( 

sm 


MTFfit(f)  = 


7lL 

fc 


x-f 

fc 


=  sme 


Kfc  J 


(4.5) 


where : 

MTFfit(f)  =  Fitting  curve  for  the  MTF 
f  =  spatial  frequency 
fc  =  cutoff  frequency 
a  =  fitting  parameter 

For  the  specific  set  of  data  collected  for  the  AGA  780 
in  the  configuration  without  polarizer,  the  best  value  for 
"a"  was  1.17.  Figure  4.4  also  shows  the  fitting  curve. 
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Figure  4.4  -  MTF  of  the  AGA  780  (no  polarizer)  Fitted  with  the 

"sine"  Function  of  Equation  4.4. 

D.  AGA  MRTD 

1.  Unpolarized  Case 

Using  Matlab®  to  calculate  and  plot  the  MRTD  from 
Equation  3.27  using  known  or  estimated  AGA  parameters  for 
the  unpolarized  case,  we  obtained  the  curve  shown  in  Figure 
4.5.  Table  4.2  shows  all  the  values  used  in  this 
calculation  with  their  respective  sources  (when  available) 
or  the  rationale  for  their  estimation. 
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Parameter 

Value 

Unit 

Source/Rationale 

SNRt 

6.5 

In  order  to  have  a 
probability  of 
detection  close  to 

100%. 

D 

0.055 

m 

Operating  Manual 

D* 

4 . 13el0 

Cm  Hz*  W'-1 

Calculated  from  the 
given  NET 

N 

1 

- 

Operating  Manual 

^Isc 

0.75 

Typical  value 

To 

0.7 

Typical  value 

a 

1.1 

mrad 

Operating  Manual 

B  1 

1.1 

mrad 

Operating  Manual 

MTF  system 

- 

— 

Measured  in  the 
laboratory 

L 

7 

— 

Standard  pattern  was 
used 

Te 

0.2 

S 

Typical  value  given 
by  Shumaker  [Ref.  8, 
page  8-38] 

dN 

dT 

6 . 27e-5 

(watt  cm-2 

K'1  sr"1) 

Shumaker  [Ref.  8, 
page  2-36] 

FOVx 

7 

Degrees 

Operating  Manual 

_ 

7 

Degrees 

Operating  Manual 

Table  4.2  -  MRTD  Parameters  with  no  Polarizer.  "Most 
Appropriate"  Set  Evaluated  for  the  AGA  780. 
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0  0.1  0.2  0.3  0.4  0.5 

spatial  frequency  in  mrad 

Figure  4.5  -  AGA  780  MRTD  (Unpolarized  case)  Based  on 
Estimated  Parameters  Shown  in  Table  4.2. 


2 .  Polarized  Case 


For 

the  polarized 

case, 

the 

same 

procedure 

was 

adopted, 

except  that  a  new  set 

of 

data 

for  the  MTF 

was 

taken. 

and  the  formula 

for 

the 

MRTD 

calculation 

was 

modified  as  shown  in  Equation  3.32.  Table  4.3  shows  the  new 
values  used  in  the  calculation  of  the  MRTD  for  the 
polarized  case,  with  their  respective  sources  (when 
available)  or  the  rationale  for  their  estimation. 


Parameter 

Value 

Unit 

Source /Rationale 

SNRt 

6.5 

In  order  to  have  a 
probability  of 
detection  close  to 

100%. 

D 

0.055 

m 

Operating  Manual 

D* 

4 . 13el0 

warn 

Calculated  from  the 
given  NET 

N 

1 

Operating  Manual 

^Isc 

0.75 

- 

Typical  value 

to 

0.7 

- 

Typical  value 

V 

0.85 

;  — . 

Polarizer 

specification 

a 

1.1 

mrad 

Operating  Manual 

3 

1.1 

mrad 

Operating  Manual 

MTFpolarizer 

•’  ^  ;  |i 

Estimated  from  the 
measured  MTF  of  the 
system  with  the 
polarizer 

L 

7 

-* 

Standard  pattern  was 
used 

Te 

0.2 

S 

Typical  value  given 
by  Shumaker  [Ref.  8, 
page  8-38] 

FOVx 

7 

Degrees 

Operating  Manual 

FOVy 

7 

Degrees 

Operating  Manual 

k 

0.5 

i'iiSS  .  ■  =1 

Blackbody 

approximation 

BN 

BT  : 

4.00e-5 

(watt  cm- 2 
K-l  sr-1) 

;|ipf  ^  | 

Matlab®  program  dNdT 
(Appendix  A) 

T 

Table  4.3  -  MRTDp  Parameters.  The  Shading  Indicates  the 
Changed  Parameters  in  Relation  to  Table  4.2. 


The  plot  of  the  MRTD  for  the  polarized  case  in 
comparison  to  the  unpolarized  configuration  is  presented  in 
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Figure  4.6.  As  expected,  the  curve  representing  the  system 
with  the  polarizer  was  shifted  upwards. 


0  0.1  0.2  0.3  0.4  0.5 

spatial  frequency  in  mrad 


Figure  4.6  -  AGA  780  MRTD  (Polarized  vs.  Unpolarized). 
Comparison  of  the  MRTDp  Computed  from  Parameters  of  Table  4.3 

with  MRTD  from  Table  4.2. 
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V.  ANALYSIS  OF  THE  RESULTS 

A.  THE  RESOLUTION  PROCESS 

While  the  preliminary  data  measurements  were  being 
taken  in  the  laboratory,  some  characteristics  of  the 
process  of  resolving  a  four-bar  pattern  for  the  MRTD 
measurements  turned  out  to  be  very  interesting.  Some  of 
them  inspired  the  study  of  the  dynamics  of  this  process  and 
the  results  derived  from  this  part  of  the  work  will  be 
briefly  discussed  in  this  section. 

The  question  to  be  answered  when  measuring  the  MRTD  is 
at  what  point  is  the  four-bar  pattern  actually  seen  as  a 
four-bar  pattern.  Since  it  is  a  subjective  decision,  some 
factors  have  great  influence  on  the  outcome: 

•  How  well  "trained”  is  the  observer? 

•  What  is  the  dynamic  of  the  experiment?  (Does  temperature 
rise  in  discrete  steps  and  it  can  be  stabilized  at  a 
precise  value?) 

•  Is  there  any  delay  between  the  time  at  which  the  observer 
identifies  the  four-bar  pattern  and  the  value  is 
recorded? 
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•  What  are  the  similarities  among  the  recorded  images  of 
the  objects  at  the  MRTD? 

The  transition  between  an  unresolved  pattern  and  a 
"barely"  resolved  one  should  be  fully  understood  in  order 
to  make  sure  the  value  being  read  is  the  actual  MRTD.  The 
general  shape  of  the  MRTD  can  be  illustrated  by  Figure  5.1 
(compare  with  Figure  4.5).  The  figure  is  divided  into  four 
basic  regions. 


Spatial  Frequency 

Figure  5.1  -  MRTD  General  Shape  and  Parameter  Space  Regions. 
Temperature  Difference  (At)  Unresolved. 
Temperature  Diference  Detectable  but  not  Resolvable. 
Temperature  Difference  Resolvable  but  Operator  Dependent. 
Temperature  Clearly  Resolved. 

Region  A  identifies  the  parameter  pairs  (spatial 

frequency  and  temperature  difference)  for  which  the  target 

(four-bar  pattern)  can  not  be  resolved,  mainly  because  of 

the  effect  of  the  Modulation  Transfer  Function  causing 

attenuation  at  the  higher  frequencies .  Region  B  represents 
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the  pairs  which  are  not  yet  resolved  by  the  observer, 
although  thorough  analysis  of  the  output  would  reveal  that 
the  four-bar  characteristic  is  present.  However,  the 
temperature  difference  between  peaks  and  valleys  (Figure 
5.2)  in  the  image  is  less  the  minimal  temperature 
difference  that  the  system  (operator  included)  can  discern. 
Region  C  represents  a  region  where  MRTD  measurements  taken 
by  different  observers  or  the  same  observer  in  different 
circumstances  may  fall.  It  is  easy  to  notice  that  it  would 
be  extremely  difficult  to  establish  a  monotonic  and  smooth 
line  to  represent  the  MRTD  by  laboratory  measurements  in 
which  the  subjectiveness  of  the  human  operator 
interpretation  plays  an  important  role.  Finally,  region  D 
corresponds  to  well  resolved  images  (low  frequency/high 
temperature  difference  pairs) . 


Peak 


Four-bar  pattern 


Equivalent  square  wave 
function  (Scan  of  pattern) 


Figure  5.2  -  Four-bar  Pattern  (Peaks  and  Valleys)  with 
"in-scan"  Xrradiance  Profile. 
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Output  examples  of  the  four  regions  can  be  viewed  in 
Figure  5.3,  which  represents  a  simulation  run  in  Matlab® 
(sim_mrt2.m.  Appendix  A).  The  region  of  most  interest  in 
this  study  is  of  course  the  boundary  between  regions  B  and 
C  (MRTD)  ,  and  especially  region  C,  where  most  of  the 
measurements  are  expected  to  be  found.  The  next  section 
will  propose  a  method  for  calibrating  values  in  the  region 
C,  so  that  they  can  better  represent  the  actual  MRTD,  and 
will  be  standardized  between  observers. 


Region  C  Region  D 


Figure  5.3  —  Region  A  Represents  an  Unresolved  Spatial 
Frequency.  In  Region  B,  the  Four-Bar  Structure  is  Present, 
but  the  Operator  can  not  Resolve  the  Bars.  Regions  C  and  D 
Show  Well  Resolved  Spatial  Frequencies. 
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B.  A  PROPOSED  METHOD  OF  CALIBRATION 

As  explained  in  the  previous  section,  region  C  is  of 
great  importance  since  many  of  the  measurements  will  be 
taken  there.  The  main  reasons  are  the  psychophysics  of  the 
human  recognition  process  and  the  dynamics  of  the 
experiment.  The  solution  proposed  to  reduce  these  effects 
and  make  use  of  the  measurements  taken  in  region  C  is  to 
establish  a  standard  output  that  would  be  considered  an 
optimum  MRTD  measurement .  Analyses  were  done  on  several 
MRTD  data  sets  and  the  conclusion  reached  was  that  the 
observer  would  detect  the  four-bar  pattern  as  soon  as  the 
temperature  difference  between  the  peaks  and  valleys  in  the 
output  was  equal  to  the  minimal  discernable  temperature  of 
the  system.  Figure  5.4  shows  the  output  of  an  ideal  MRTD 
measurement.  Notice  that  the  temperature  difference  between 
valleys  and  peaks  is  approximately  0.1°C,  which  is  the 
minimal  discernable  temperature  for  the  AGA  780  [Ref.  4]. 
It  is  important  to  notice  that  the  usual  situation  is  of 
Tt>Tb,  hot  target  against  cool  background  -  or  "white  hot." 
The  situation  chosen  in  this  analysis  (Tb<Tt)  represents  a 
reversed  contrast,  and  was  chosen  due  to  the  simplicity  of 
representing  it  in  the  laboratory. 
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Figure  5.4  -  Ideal  MRTD  Measurements.  Detected  Signal  vs. 

Angular  Position,  Scan  Direction. 

In  Figure  5.4,  we  have : 

Tt  =  target  temperature 

Tb  =  background  temperature 

Tvaiiey  =  apparent  temperature  in  the  valley 

Tpeak  =  apparent  temperature  of  the  peak 

Based  on  this  "ideal"  output  and  the  observation  of 
the  typical  output  falling  in  region  C,  we  can  establish  a 
relationship  between  those  values.  By  doing  so,  we 
"calibrate"  the  MRTD  measurements  in  region  C  to  more 
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accurate  values.  This  calibration  is  illustrated  in  Figure 


5.5. 


Figure  5.5  -  Calibration  Process  for  Observed  MRTD.  The  Value 
of  tmeasurement-tMRTD  Corresponds  to  the  Time  Lag  between  Ideal 

and  Non- Ideal  Measurements. 

The  experiment  was  done  so  that  the  target  temperature 
was  held  constant  and  the  background  temperature  was 
increased,  and  the  observer  would  determine  when  the  four- 
bar  pattern  would  be  just  resolved.  A  delay  in  determining 
the  actual  instant  when  the  MRTD  was  reached  would  lead  to 
a  measurement  falling  in  region  C,  which  has  a  greater 
value  of  Tb  and  Tvaixey.  In  addition,  the  difference  between 
Tb  and  Tvaiiey  is  greater  than  the  minimum  discernable 
temperature  (Tmd)  . 
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Using  properties  of  similar  triangles  we  can  derive  a 
formula  for  Tc  (calibrated  temperature) .  Since  triangles  ABC 
and  ADE  are  similar. 


DE  =  DG^  Tb  -rv  =  7 1 -Tt 
BC  BF^  ATmd  Tc-Tt 


from  which 

AT.,-(r»-r,)+r,-(rt-r,) 

‘  (Tt-Tr) 


so  that 


MRTDcal 


=  T  -T  = 

x  C  Xt 


AT^-CTt-T,) 


Tb~Tv 


(5.3) 


where : 

ATmd  =  minimum  discernable  temperature 
Tc  =  calibrated  background  temperature 
MRTDcai  =  calibrated  MRTD. 

Figure  5.6  shows  a  comparison  between  a  data  set  taken 
by  an  untrained  observer  and  the  corresponding  calibrated 
version.  As  a  reference,  the  calculated  value  of  the  MRTD 
is  also  included.  The  uncalibrated  data  show  considerable 
scatter,  with  a  "fitted"  value  that  deviates  widely  from 
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the  calculated  MRTD.  The  "calibrated"  data  however  show 


much  reduced  scatter  and  a  fitted  value  in  much  closer 
agreement  with  the  calculation.  For  the  trained  observer 
the  measured  and  calibrated  data  show  very  close  agreement . 
This  means  that  the  trained  observer  is  much  more  likely  to 
determine  the  exact  moment  when  the  four-bar  pattern  is 
just  resolvable.  Figure  5.7  shows  an  equivalent  data  set  to 
that  of  5.6,  but  taken  by  a  trained  observer.  Although 
Figures  5.6  and  5.7  also  represent  different  polarization 
filters,  it  will  not  affect  the  analysis  of  the  improvement 
in  quality  of  the  data  collected  from  a  trained  observer. 
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Figure  5,6  -  Comparison  of  Uncalibrated,  Calibrated  and 
Calculated  MRTD  (Untrained  Observer  -  Unpolarized  Case) . 
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Figure  5-7  -  Comparison  of  Uncalibrated  and  Calibrated  MRTD 
for  a  Trained  Observer  (Vertical  Polarization  Filter) . 

C.  CALCULATED  DATA  VERSUS  MEASURED  DATA 

After  collecting  the  laboratory  data,  the  next  step 
was  to  verify  how  close  the  measured  values  were  to  the 
calculated  ones.  Using  the  MRTD  developed  in  Chapter  III 
and  the  parameters  in  Tables  4.2  and  4.3,  we  compared 
calculated  and  measured  data.  Figure  5.8  shows  the 
comparison  for  the  unpolarized  case.  In  this  case,  three 
sets  of  data  were  taken  in  different  days  (May  5th,  11th,  and 
25th)  .  We  can  observe  the  scatter  decreasing  with  repetition 
of  the  task,  with  the  third  set  (May  25th)  corresponding  to 
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the  trained  observer,  showing  very  little  scatter  and 
excellent  agreement  with  the  calculation. 


0.1  0.2  0.3  0.4 

spatial  frequency  in  cycles/m  rad 


0.5 


Figure  5.8  -  Comparison  of  Measured  and  Calculated  MRTD 
(Unpolarized  Case),  and  Influence  of  Training. 

X:  Measurement  of  May  5; 

:  Measurement  of  May  11; 

O:  Measurement  of  May  25. 

Figures  5 . 9  and  5 . 10  show  equivalent  comparisons  to 
data  of  Figure  5.8  but  taken  using  polarization  filters. 
Filters  4  and  5  correspond  to  the  horizontal  and  vertical 
polarization  filters  respectively.  In  both  plots,  the 
comparison  is  made  between  "trained  observer"  data  and  the 
calculated  MRTDP. 
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spatial  frequency  in  cycles/m  rad 


Figure  5.9  -  Comparison  of  Measured  and  Calculated  MRTDp 
(Horizontal  Polarization  Case  -  Filter  4) . 


spatial  frequency  in  cycles/m  rad 


Figure  5.10  -  Comparison  of  Measured  and  Calculated  MRTDp 
(Vertical  Polarization  Case  -  Filter  5) . 
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D.  RANGE  CALCULATION 

The  final  analysis  in  this  thesis  is  the  range 

calculation.  The  maximum  range  at  which  some  specific  task 
can  be  performed  is  estimated  using  the  apparent  target- 
background  temperature  difference  (TBTD)  and  the  MRTD 
curves,  developed  in  Chapter  III.  Although  MRTD  is  a 
function  of  spatial  frequency,  we  can  translate  it  into  a 
function  of  range,  if  specific  target  and  task  are  given. 
Some  criteria  have  been  established  to  relate  the 
resolution  of  bar  charts  to  levels  of  visualization  [Ref. 
8].  Johnson  [Ref.  10]  first  suggested  his  criteria  in  1957. 

Others  such  as  Moser  and  O'Neill  followed  him  with  their 

own  criteria.  Basically,  all  the  criteria  try  to  establish 
the  number  of  resolution  elements  (resels)  or  line  pairs 

per  critical  dimension  needed  to  perform  a  given  task.  It 
is  then  possible  to  convert  spatial  frequency  into 
distance,  as  seen  in  Figure  5.11.  It  is  important  to  notice 
that  one  spatial  cycle  corresponds  to  two  resels .  We  need 
then  to  have  a  full  cycle  of  spatial  frequency  on  the 
average  target  dimension  being  observed  to  achieve 
detection.  This  dimension  is  also  known  as  "critical 
dimension,"  and  is  defined  differently  by  various  authors. 
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Dc 


Figure  5.11  -  Conversion  of  Spatial  Frequency  into  Distance 

(Task:  Detection) . 

Using  the  small  angle  approximation,  and  applying  the 
appropriate  factor  to  get  the  result  in  cycles  per 
milliradians,  we  have  the  following  formula: 

R  _  2000 DJ_  (5.4) 

N 

where : 

R  =  Range  (m) . 

Dc  =  Critical  Dimension  (m) . 
f  =  spatial  frequency  ( cycles /mrad) . 

N  =  Number  of  resels  per  critical  dimension  (Dc)  ,  based 
on  some  pre-established  criteria.  For  example,  N  would  be 
equal  to  two  for  detection. 
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Figure  5 . 12  shows  the  unpolarized  and  polarized  MRTD 
developed  in  Chapter  III  plotted  as  a  function  of  range  for 
detection  (two  resolution  elements  in  the  critical 
dimension)  of  a  target  with  critical  dimension  of  15m 
(comparable  to  the  Research  Vessel  "POINT  SUR"  used  in  the 
EOPACE3  campaign) . 


Range  (km) 

Figure  5.12  -  MRTD  for  the  AGA  780  as  a  Function  of  Range 
(Polarized  and  Unpolarized  Cases) 


3  EOPACE  -  Electro  Optic  Propagation  Assessment  in  the  Coastal  Environment  was  a 
multinational  measurement  series  in  San  Diego  Bay  in  1996,  sponsored  by  ONR,  and 
organized  by  SPAWAR  SYSTEMS  CENTER.  Its  overall  purpose  was  to  quantify  infrared  (IR) 
propagation  characteristics  for  near  ocean  surface  transmission  and  analyze  electro-optic 
(EO)  systems  performance  in  the  coastal  environment  condition  [Ref.  12] 
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Now,  it  is  necessary  to  model  the  second  part  of  the 
problem,  the  apparent  temperature  difference  as  a  function 
of  range.  The  main  interest  of  this  thesis  is  in  the 
maritime  environment,  where  the  sea  surface  emission 
(background)  near  the  Brewster  angle  in  the  8-12(im  region 
has  been  shown  to  exhibit  strong  vertical  polarization 
[Ref.  1].  By  using  a  horizontal  polarizer,  it  is  possible 
to  filter  out  some  of  the  background  radiation,  thus 
improving  the  contrast. 

The  observed  degree  of  polarization  of  the  radiance  is 
defined  as  [Ref.  1]: 


POL  = 


(».)+(»*) 


(5.5) 


where : 

<NV>  =  mean  apparent  vertically  polarized  radiance. 

<Nh>  =  mean  apparent  horizontally  polarized  radiance. 

It  means  that  a  degree  of  polarization  greater  than 
zero  implies  a  net  vertical  polarization,  while  zero 
corresponds  to  an  unpolarized  source.  Table  5.1  tabulates 
the  percentages  of  vertical  and  horizontal  polarization 
corresponding  to  some  values  of  POL.  The  mathematical 
analysis  involved  in  this  table  is  carried  out  in 
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Equations  5.6  to  5.12.  We  define 


W  ma  W 

<JV.)  +  (JV()  (JV.)+(JV,> 


(5.6) 


as  the  fractions  of  vertical  and  horizontal  polarization 
radiance.  If  we  subtract  those  expressions,  we  have: 


{N,)+{Nt)  (N.)+(N„)  {N.)+{N„) 


(5.7) 


which  is  the  same  as  equation  5.5.  Adding  the  two 
expressions  in  5.6  the  result  is  one. 


W  +  W  W+M  x 

W+K)  (*,)+<*»}  (JV„>  +  <\,) 


(5.8) 


We  then  have  a  system  with  two  equations  and  two 


unknowns : 


=  POL 


( N,)+(Nh )  '  (N,)+(N„) 

(«■)  ,  W  ;1 

(*,)+<*.)  (w,)+<w,) 


(5.9) 


(5.10) 


Solving  the  system  above  for  the  fractions  of  vertical 
and  horizontal  polarization  radiance,  we  have 


K) 

<N»>  +  K> 

fa.) 

(Hv>  +  <Nh) 


(l  +  POL) 


(l  -  POL) 


:5.11) 


(5.12) 
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POL 

Percent  Vertical 
polarization 

Percent  Horizontal 
polarization 

0  1 

50% 

50% 

0.1 

55% 

45% 

0.2 

60% 

40% 

0.3 

65% 

35% 

POL  H 

( 1+POL) /2 

(1-POL) /2 

Table  5.1  -  Degrees  of  Polarization  and  Respective 

Percentages . 


Typical  values  for  the  degree  of  sea  polarization 
measured  during  the  MAPTIP4  campaign  range  from  7  to  30%  in 
the  LWIR  band.  Assuming  the  target  is  unpolarized  (POL=0), 
we  can  use  a  horizontally  polarized  filter  to  eliminate 
some  of  the  background  radiation.  Knowing  target  and 
background  temperatures,  we  can  calculate  the  mean  radiance 
at  zero  range  (Equations  2.5  and  3.2)  .  A  perfect 
horizontally  polarizing  filter  would  eliminate  50%  of  the 
target  radiance  and  (l+POL)/2  of  the  background  radiance 
(see  Table  5.1).  Since  only  a  fraction  of  the  total 
radiance  is  passing,  the  sensor  "sees”  the  scene  as  if  the 
target  and  background  both  had  lower  temperatures.  It  means 
that  the  observed  temperature  is  smaller  than  the 
emitted.  The  observed  temperature  can  be  calculated  by 
solving  Equation  3.28  for  TB,  when  "k"  is  equals  to  the 


4  MAPTIP  -  Marine  Aerosol  Properties  and  Thermal  Imager  Performance  Trial  campaign  was 
held  in  coastal  waters  of  Katwijk,  Netherlands,  in  1993.  Its  objective  was  "to  improve 
marine  aerosol  models  for  the  marine  boundary  layer  and  the  modeling  of  electromagnetic 
propagation  and  imaging  in  the  coastal  environment."  [Ref.  1] 
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fraction  of  the  total  radiance  passing  through  the  filter, 
in  this  case  (l-P0L)/2.  An  apparent  target-background 
temperature  difference  (TBTD)  at  zero  range  can  then  be 
established. 


Table  5 . 2  shows  some  values  of  apparent  TBTD  for  the 
case  where  target  temperature  (Tt)  is  310°  and  background 
temperature  (T.B)  is  300°K.  Appendix  D  shows  a  table  of 
values  for  <N>  for  temperatures  ranging  from  240  to  320°K. 


POL 

7% 

10% 

20% 

30% 

<Nt>  (310°) 

1 . 1278e-3 

1 . 1278e-3 

1 . 1278e-3 

1 . 1278e-3 

<Nb>  (300°) 

9 . 6243e-4 

9 . 6243e-4 

9 . 6243e-4 

9 . 6243e-4 

<Nt'>  with 

polarizer 

5 . 639e-4 

5 . 639e-4 

5 . 639e-4 

5 . 639e-4 

<Nb'>  with 

polarizer 

4 . 4753e-4 

4 . 3309e-4 

3 . 8497e-4 

3 . 3685e-4 

Tt' 

270. 4°K 

270. 4°K 

270. 4°K 

270. 4°K 

Tb' 

259. 2°K 

257. 7°K 

252. 4°K 

246. 7°K 

AT  (unp.) 

O 

o 

H 

H1 

O 

o 

o 

0 

o 

0 

AT'  (pol.) 

11.2° 

12.7° 

o 

00 

r-H 

23.7° 

Table  5.2  -  Calculation  of  Apparent  TBTD  as  a  Function  of  the 
Degree  of  Polarization  of  the  Background  with  Unpolarized 

Target . 

In  Table  5.2,  the  following  parameters  are  shown: 

<Nt>  =  mean  emitted  target  radiance. 

<Nb>  =  mean  emitted  background  radiance. 

<NT>  =  mean  observed  target  radiance  (<NT>*0.5). 

<Nb  >  =  mean  observed  background  radiance  ( (<NT>*  (1-POL)  /2)  . 

Tt '  =  apparent  target  temperature  (from  Equation  3.28). 
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Tb'  =  apparent  background  temperature  (from  Equation  3.28)  . 
AT  =  actual  temperature  difference  (310°-300°K) . 

AT'  =  observed  temperature  difference  (Tt'-Tb'). 

The  next  step  is  to  calculate  the  apparent  temperature 
as  a  function  of  range,  using  a  computer  code  that 
estimates  atmospheric  transmittance,  such  as  LOWTRAN, 
MODTRAN,  or  SEARAD.  A  Mat lab  program  (rangeimp.m  -  Appendix 
A)  was  written  to  calculate  and  plot  the  MRTD  and  the 
apparent  TBTD  versus  range  for  the  example  presented  in 
this  section.  This  program  takes  the  apparent  temperature 
difference  shown  in  Table  5.2  and  multiplies  it  by  a 
transmissivity  factor  obtained  from  SEARAD  for  the 
following  set  of  input  parameters: 

•  Lowtran7  was  selected. 

•  Type  of  atmospheric  path:  Slant  path  between  two 
altitudes . 

•  Surface  albedo  of  earth:  assumed  blackbody. 

•  Navy  maritime  model. 

•  No  clouds  or  rain. 

•  Altitude  of  the  target:  0. 

•  Altitude  of  the  sensor:  300ft. 

•  Wind  speed:  7 . 5  m/s . 
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•  Wind  speed:  7.5  m/s  (24h  average). 


•  Airmass  character:  3. 

•  Visibility:  the  program  calculated  visibility  from  the 
standard  data  for  the  chosen  conditions  (Mid-latitude 
summer) . 

Figures  5.13  and  5. 14  (zoom)  show  the  result  of  this 
calculation,  which  is  summarized  in  Table  5.3.  Figure  5.13 
shows  the  apparent  TBTD  computed  for  various  degrees  of 
polarization  as  a  function  of  range  for  the  sample 
atmospheric  conditions.  The  MRTD  for  detection  with  and 
without  the  polarizer  is  plotted  on  the  same  scale  for  an 
assumed  15m  target  critical  dimension.  For  any  case  the 
intersection  of  corresponding  TBTD  and  MRTD  curves  defines 
the  detection  range.  Figure  5.14  shows  a  magnified  plot  of 
the  intersection  region  of  Figure  5.13. 

We  can  observe  that  for  degrees  of  polarization 
greater  then  20%  we  obtain  an  improvement  in  range  by  using 
the  polarizing  filter.  On  the  other  hand,  values  of  POL 
smaller  than  20%  imply  degraded  range  performance,  since 
the  effects  of  the  polarizer  in  increasing  MRTD  are  more 
significant  than  the  improvement  in  the  apparent  TBTD. 
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Range  (km) 

Figure  5.13  -  MRTD  and  TBTD  as  Functions  of  Range,  for  Various 
Degrees  of  Polarization.  MRTD  is  for  Detection  with  Critical 
Dimension  15m.  Intersection  Points  Define  Maximum  Ranges  for 

Detection. 


Figure  5.14  -  MRTD  and  TBTD  as  Functions  of  Range  ( Zoom)  . 


K 


POL 

Range  (km) 

Difference  (km) 

0% 

11.13 

0 

10% 

10.25 

00 

00 

• 

o 

1 

20% 

11.06 

1 

o 

• 

o 

-4 

30% 

11.63 

0.60 

Table  5.3  -  Difference  in  Range  Performance  for  Increasing 

Degrees  of  Polarization. 


All  the  calculation  done  so  far  approximates  the  real 
world,  since  it  assumes  a  simple  model  where  only  the 
radiance  from  target  and  background  are  accounted  for.  More 
elaborate  models  have  been  proposed,  such  as  the  one  shown 
in  Lagaras  [Ref.  11],  which  uses  experimental  data  from  the 
EOPACE  campaign.  In  his  work,  Lagaras  suggested  a  TBTD  for 
both  eases,  with  and  without  polarization,  using  target 
temperature  of  300.26°K  and  background  temperature  of 
289.95°K,  as  shown  in  Figure  5.15. 
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Range  (km) 


Figure  5.15  -  TBTD  for  Polarized  and  Unpolarized  Cases  [After 

Ref.  11] . 

.Comparing  this  new  set  of  TBTD  curves  and  our 
formulation  for  the  MRTD  for  both  cases,  we  notice  a  small 


improvement  in  range,  as  shown  in  Figure  5.16. 


25 


Figure  5.16  -  MRTD  and  TBTD  as  Functions  of  Range. 


Both  results  presented  in  Figures  5.14  and  5.16  show 
that  it  is  possible  to  achieve  some  improvement  in  range 
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performance  for  thermal  imaging  systems.  However  the 
results  also  pointed  out  to  a  relatively  small  improvement, 
if  any.  It  should  be  noted  that  the  calculations  were  made 
for  only  a  single  target  critical  dimension,  a  single 
atmosphere  state,  and  two  TBTD  cases. 
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VI .  CONCLUSION 


The  main  objective  of  this  thesis  was  to  determine  a 
more  appropriate  Minimum  Resolvable  Temperature  Difference 
(MRTD)  formulation  for  use  with  a  polarization  filtered 
thermal  imaging  system  (generic  FLIR) .  In  order  to  analyze 
the  problem,  two  approaches  were  used. 

First,  a  mathematical  formulation  was  derived  from  the 
traditional  MRTD  expressions  used  by  various  authors  [Ref. 
8  and  Ref.  9].  Polarization  filter  transmissivity, 
Modulation  Transfer  Function  (MTF) ,  and  the  "thermal 
derivative"  of  Planck's  Equation  were  studied  and  their 
effects  introduced  to  the  new  MRTD  formulation,  when 
applicable.  In  a  second  phase,  laboratory  experiments  using 
the  AGA  Thermovision  780  system  were  conducted.  The  MRTD 
patterns  were  set  up  (see  Appendix  B)  with  the  object  of 
gathering  data  to  be  compared  to  the  proposed  mathematical 
formulation. 

The  main  conclusions  of  this  work  are  listed  below: 

•  The  mathematical  formulation  suggests  that  use  of 
polarization  filters  increases  the  temperature 
difference  at  which  the  observer  detects  the  four- 
bar  pattern.  This  reduction  in  performance  depends 
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on  the  polarization  filter  characteristics  as  well 
as  on  the  new  reference  temperature  to  be  used  in 
the  "thermal  derivative"  of  Planck's  Equation. 

•  The  experimental  data  agreed  with  the  mathematical 
formulation.  Data  collected  in  the  laboratory  showed 
very  small  deviation  from  the  proposed  MRTD 
developed  in  this  thesis.  Figure  6.1  summarizes  both 
the  results  of  the  mathematical  formulation  and  the 
laboratory  results . 


spatial  frequency  in  cycles/mrad 

Figure  6.1  -  Summary  of  Mathematical  and  Experimental  Values 
of  MRTD  with  and  without  Polarizing  Filters. 
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•  The  process  of  resolving  the  four-bar  pattern  (or 
any  type  of  target)  is  very  complex.  Training  plays 
a  decisive  role  in  the  ability  to  determine  the 
exact  MRTD  value.  In  Chapter  IV,  a  calibration 
method  was  proposed  to  adjust  measurements  taken  by 
an  untrained  observer  or  under  non-ideal  conditions. 
The  use  of  this  method  would  enable  the  untrained 
observer  to  attain  results  comparable  to  the  better- 
trained  one . 

•  Range  calculation  suggested  that  the  use  of 

polarization  filters  might  improve  the  performance 
of  thermal  imaging  systems  in  a  marine  environment, 
under  certain  conditions.  The  main  factor  to  be 
determined  is  the  degree  of  polarization  of  the  sea 
surface.  A  high  degree  of  polarization  would  lead  to 
an  improvement  in  range.  On  the  other  hand,  low 
degrees  of  polarization  would  imply  degradation  of 
the  system  performance.  The  example  used  in  this 
research  showed  improvement  in  range  when  the  degree 
of  polarization  was  greater  than  20%. 

It  is  important  to  emphasize  that  this  work  is  not 
self-contained.  Extensive  work  was  previously  done  in 
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collecting  field  data  as  well  as  modeling  the  target  and 
the  background  [Ref.  1,  11,  12,  13,  and  14].  There  are 

still  some  aspects  to  be  analyzed  or  improved.  As  a 
suggestion  for  future  work,  it  would  be  interesting  to 

refine  the  target- to-background  temperature  models.  It 
would  allow  better  range  improvement  calculations.  In 
addition,  more  precise  laboratory  measurements  would  help 
to  validate  the  mathematical  formulation  proposed  here. 
This  would  be  necessary  to  have  a  temperature-controlled 
laboratory,  with  precise  instruments,  and  an  adequate 
number  of  trained  observers. 

The  results  of  this  work,  added  to  the  results  of 
others  in  this  area,  may  be  useful  in  improving  the 

performance  of  thermal  imaging  systems.  The  expected 

improvement  margin,  though  not  large,  may  be  enough  to  help 
systems  with  limited  range  capability.  Moreover,  tactical 
decision  assessment  for  electro-optic  systems  operating  in 
a  marine  environment  should  include  results  of  this  work. 
It  would  allow  them  to  achieve  higher  effectiveness  when 
assessing  decisions  over  situations  where  the  maximum 
possible  performance  is  expected. 
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APPENDIX  A.  MATLAB  SOURCE  CODES 


The  following  Matlab  codes  are  presented  in  this 
Appendix : 


FILE  NAME 

CHAPTER  OF 

REFERENCE 

f ig2_2 .m 

II 

f ig2_3 .m 

II 

f ig4_3 .m 

IV 

mtf f it .m 

IV 

mrtd_plot .m 

IV 

sim_mrt2 .m 

V 

dl2temp .m 

V 

plotcal .m 

V 

plotcal5.m 

V 

f ig5_8 .m 

V 

f ig5_9 .m 

V 

f ig5_10 .m 

V 

dndt . m  ; 

V 

Enne .  m 

V 

Rangeimp4 ,m 

V 

Rangeimp5 .m 

V 

Fig6_l .m 

VI 

readptw2 . m 

IV5 

readmtf .m 

IVs 

f  ilt  .m 

IV7 

guil .m 

Tv8 

readgui .m 

IV9 

Table  A.l  -  List  of  Matlab®  Programs. 


5  ,6,  7,  B,  and  9  - 


those  programs  were  not  referenced  in  Chapter  IV,  but  were  used  to  analyze 
most  of  the  results  collected  during  the  experiments  discussed  in  that  chapter. 
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<^********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  fig2_2 

%  Function:  plots  Planck's  Equation  for  four  different 
%  temperatures  and  shows  Wien's  Law 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  calculates  Planck's  Equation  for  T=220,  250,  280, 

%  300,  320°K 

%  2 )  plots  the  results 

%********************************************************** 

^**j'*************************‘k*mk'k**',k**'k**'k'k'k*'k-k-k'k-*:'k'k'k'k'k-k'k-k'*: 


clear 
close  all 
lambda=0 : 0 . 01 : 25  ; 

%  Planck's  Equation  for  T=220 
T0=22  0 ; 
cl=3 . 7418e4 ; 
c2  =  l . 43  88e4 ; 

M0=exp (c2 . / (lambda*T0) ) -1; 

Ml= (cl . / (lambda. ^5) ) ; 

M2  2  0  =M1 . /M0 ; 

%  Planck's  Equation  for  T=250 
Tl=250 ; 

M2=exp (c2 . / (lambda*Tl) ) -1; 

M2  5  0  =Ml . /M2 ; 

%  Planck's  Equation  for  T=280 
T2=280 ; 

M3=exp (c2 . / (lambda*T2 ) ) -1; 
M280=M1. /M3; 

%  Planck's  Equation  for  T=300 
T3=3  00 ; 

M4=exp (c2 . / (lambda*T3) ) -1; 

M3  00=M1 . /M4 ; 

%  Planck's  Equation  for  T=320 
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T4=320; 

M5=exp(c2 . / (lambda*T4) ) -1; 

M3  2  0  =M1 . /M5 ; 

%  Wien's  Displacement  Law 
10=2897/T0; 

11=2897/T1; 

12=2897 /T2 ; 

13=2897/T3 ; 

14=2897/T4; 

L=  [10, 11,12,13, 14] 

maxima=  [max(M220)  ,max(M250)  ,max(M280)  ,max(M300)  ,max(M320)  ]  ; 

%  plots  the  results 
figure (1) 

plot ( lambda, M250, lambda, M2 80, lambda, M3 00, lambda, M3 20 ,  L, 

maxima, ' : ' ) 

grid  on 

gtext ( ' 250  K' ) 

gtext ( ' 280  K' ) 

gtext ( ' 300  K' ) 

gtext ( '320  K' ) 

xlabel ( 'Wavelength  -  micrometer') 
ylabel (' Spectral  Radiant  Exitance  -  M' ) 


%*********  Naval  Postgraduate  School  -  CA  **************** 


a********* 


Thesis  Research 


Type  :  Procedure 

Name  :  fig2_3 

Function:  plots  transmissivity  data  from  Searad 
Date  :  May  /  1999 

Version  :  1.0 

Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt . 
Routine  Information: 

1)  plots  the  data  collected  from  Searad 


(BAF) 


%  Midlatitude  summer 
%  Horizontal  path  lKm 
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%  Navy  Maritime  Aerosol 
%  Airmass  Characteristic  =  3 
%  No  clouds,  no  rain 
%  Sea  Level 

%  Temperature  of  the  boundary  =  30 OK 

%  lambda  =  wavelength  vector  from  2  to  20pm 

lambda= [20 . 000  19.231  18.519  17.857  17.241  16.667  16.129 
15.625  15.152  14.706  14.286  13.889  13.514  13.158  12.821 

12.500  12.195  11.905  11.628  11.364  11.111  10.870  10.638 

10.417  10.204  10.000  9.804  9.615  9.434  9.259  9.091  8.929 
8.772  8.621  8.475  8.333  8.197  8.065  7.937  7.813  7.692  7.576 

7.463  7.353  7.246  7.143  7.042  6.944  6.849  6.757  6.667  6.579 

6.494  6.410  6.329  6.250  6.173  6.098  6.024  5.952  5.882  5.814 

5.747  5.682  5.618  5.556  5.495  5.435  5.376  5.319  5.263  5.208 

5.155  5.102  5.051  5.000  4.950  4.902  4.854  4.808  4.762  4.717 

4.673  4.630  4.587  4.545  4.505  4.464  4.425  4.386  4.348  4.310 

4.274  4.237  4.202  4.167  4.132  4.098  4.065  4.032  4.000  3.968 

3.937  3.906  3.876  3.846  3.817  3.788  3.759  3.731  3.704  3.676 

3.650  3.623  3.597  3.571  3.546  3.521  3.497  3.472  3.448  3.425 

3.401  3.378  3.356  3.333  3.311  3.289  3.268  3.247  3.226  3.205 

3.185  3.165  3.145  3.125  3.106  3.086  3.067  3.049  3.030  3.012 

2.994  2.976  2.959  2.941  2.924  2.907  2.890  2.874  2.857  2.841 

2.825  2.809  2.793  2.778  2.762  2.747  2.732  2.717  2.703  2.688 

2.674  2.660  2.646  2.632  2.618  2.604  2.591  2.577  2.564  2.551 

2.538  2.525  2.513  2.500  2.488  2.475  2.463  2.451  2.439  2.427 

2.415  2.404  2.392  2.381  2.370  2.358  2.347  2.336  2.326  2.315 

2.304  2.294  2.283  2.273  2.262  2.252  2.242  2.232  2.222  2.212 

2.203  2.193  2.183  2.174  2.165  2.155  2.146  2.137  2.128  2.119 

2.110  2.101  2.092  2.083  2.075  2.066  2.058  2.049  2.041  2.033 

2.024  2.016  2.008  2.000]; 

%  t  =  transmissivity  vector  for  lambda 

t= [0 . 0025  0.0045  0.0180  0.0310  0.0310  0.0504  0.0342  0.0027 
0.0000  0.0000  0.0103  0.1287  0.2467  0.4353  0.5171  0.5290 

0.6542  0.6792  0.7016  0.7259  0.7611  0.7658  0.7870  0.7856 

0.7985  0.8166  0.7890  0.7828  0.7765  0.8055  0.7911  Q.7857 

0.7615  0.7369  0.6692  0.6575  0.5460  0.4979  0.2434  0.1874 

0.1375  0.0256  0.0083  0.0006  0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0003  0.0029  0.0054  0.0073 

0.0365  0.0912  0.2170  0.2216  0.2948  0.4341  0.4646  0.5111 

0.6017  0.6619  0.6706  0.7095  0.6769  0.4922  0.3929  0.2232 
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0.0921  0.0249  0.0001  0.0000 
0.7962  0.8291  0.8592  0.8782 
0.8365  0.8390  0.8202  0.7925 
0.7986  0.7778  0.7328  0.7102 
0.7662  0.7536  0.6786  0.5842 
0.2384  0.2441  0.1095  0.0745 
0.1270  0.1496  0.1740  0.1286 
0.2103  0.1939  0.1510  0.0780 
0.0001  0.0000  0.0000  0.0000 
0.0000  0.0000  0.0000  0.0000 
0.0000  0.0000  0.0000  0.0000 
0.1256  0.2859  0.4405  0.4001 
0.6640  0.6633  0.7131  0.7682 
0.8596  0.8581  0.8802  0.8827 
0.8908  0.8844  0.8477  0.8566 
0.8919  0.9017  0.8840  0.8620 
0.7984  0.7479  0.7332  0.7881 
0.5325  0.5743  ] ; 


0.0000  0.0000  0.0002  0.7085 
0.8972  0.8992  0.8795  0.8603 
0.7776  0.7961  0.8043  0.6836 
0.7125  0.7548  0.7826  0.7804 
0.4220  0.3888  0.2501  0.3650 
0.4278  0.5665  0.3040  0.2678 
0.2060  0.2595  0.3180  0.1514 
0.0620  0.0228  0.0036  0.0007 
0.0000  0.0000  0.0000  0.0000 
0.0000  0.0000  0.0000  0.0000 
0.0001  0.0037  0.0357  0.1192 
0.5262  0.5733  0.5078  0.6329 
0.7613  0.8134  0.8317  0.8473 
0.8818  0.8913  0.8935  0.8982 
0.8429  0.8672  0.8609  0.8502 
0.8741  0.8509  0.8541  0.8371 
0.8184  0.8024  0.7187  0.5597 


subplot (6,1, [1,2, 3, 4, 5]) 

plot (lambda,  t) 

grid 

xlabel ( 'wavelength  -  micrometer') 
ylabel ( ' Transmissivity' ) 
axis ( [2,20,0,1] ) 


9*  •k'k'k'k'k'k'k'kir 
o 


Naval  Postgraduate  School 
Thesis  Research 


%  Type  :  Procedure 

%  Name  :  fourier2 

%  Function:  compares  input  and  output  for  a  given  spatial 
%  frequency 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information : 

%  1)  inputs  the  spatial  frequency 

%  2)  calculates  the  pattern 

%  3)  process  the  input  (MTF) 

%  4 )  plots  input  and  output 

<^*****ir±**jc**1i;********1'*********i'±‘k-*:****-k'*:-k'*:‘k*-k***'k'*?*'k'k'*r'k** 
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clear 
close  all 
points=1000 ; 
mn=-50 ; 
mx=50; 

mrad=linspace (mn,mx, points)  ; 

%  gets  MTF  from  MTFfit  function 

MTFfit; 

close  all 

%  builds  a  two-side  MTF 
MTF2= [f liplr (MTFl)  MTFl] ; 

MTF3=zeros (1, length (mrad) ) ; 

MTF 3 ( (length (mrad) /2-length (MTF2 ) /2 

+1) : ( length (mrad) /2+length(MTF2) /2) )=MTF2; 

%  inputs  data 

f=input ( 'Enter  the  spatial  frequency: ' ) ; 
cycle=l/f ; 

%  builds  the  pattern 
patl=zeros(l,length(mrad)); 
uni t=cycle*points / (mx-mn) ; 

L=3 . 5*unit ; 
warning  off 

s=-cos (mrad*2*pi/ cycle) ; 
s (find(s>=0) ) =1; 
s (1 : (length (mrad) -L) /2) =0; 

s ( ( ( (length (mrad) -L) /2)+L+l) : length (mrad) ) =0; 
s (f ind(s<0 ) ) =0 ; 
warning  on 

%  FFT  of  the  pattern 
S=fft (s)  ; 

freql=0 . 5* length (mrad) / (mx-mn) ; 
freq=linspace (-freql, freql, length (mrad) ) ; 
figure (1) 
subplot (2,1,2) 

S=f ftshift (S) ; 

%  "convolution" 

OUT=MTF3 . *S; 
out=ifft (OUT) ; 
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%  plots  the  results 

figured) 

subplot (1, 2,1) 

plot (mrad, s) 

grid 

axis ( [-2* cycle  2*cycle  0  1.1]) 
subplot (1,2,2) 
plot (mrad, abs (out) ) 
grid 

axis ( [-2*cycle  2*cycle  0  1.1]) 


%********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  mtffit 

%  Function:  fits  a  curve  to  the  MTF  data 
%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  plots  the  data  collected  from  the  AGA 

%  2)  fits  the  data  to  a  'sine'  function 

%********************************************************** 


%clear 

%  MTF  measured  data 


mean0= [66 . 5592 
meanl= [66 . 5809 
mean2= [65 . 1908 
mean3= [65 . 7399 
mean4= [64.9921 
mean5= [65 . 5592 
mean6= [65.5520 


64.5289] ; 
64.5434]  ; 
63.5347]  ; 
64.7775]  ; 
64.2948]  ; 
65.5520]  ; 
65.5520]  ; 


%  calculates  modulation 

f = [ 0  0.1  0.2  0.3  0.4  0.5  0.6]  ; 

m0= (meanO ( 1 ) -meanO ( 2 ) ) / (meanO ( 1 ) +mean0 ( 2 ) ) ; 

ml= (meanl ( 1 ) -meanl ( 2 ) ) / (meanl ( 1 ) +meanl ( 2 ) ) ; 

m2= (mean2 ( 1 ) -mean2 ( 2 ) ) / (mean2 ( 1 ) +mean2 ( 2 ) ) ; 

m3= (mean3 ( 1 ) -mean3 ( 2 ) ) / (mean3 ( 1 ) +mean3 ( 2 ) ) ; 

m4= (mean4 ( 1 ) -mean4 ( 2 ) ) / (mean4 ( 1 ) +mean4 ( 2 ) ) ; 
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m5=  (mean5  (1)  -xnean5  (2)  )  /  (mean5  ( 1 )  +mean5  (2)  )  ; 
m6= (mean6 ( 1 ) -mean 6 ( 2 ) ) / (mean6 ( 1 ) +mean6 ( 2 ) ) ; 
m=  [mO  ml  m2  m3  m4  m5  m6]  ; 

MTF=m/max(m)  ; 

%  fitting  process 
min=lelO ; 

ER=  [  ]  ; 

for  x=l : 0 . 01 : 2 

MTF_fit= (sine (f / ( . 55) ) ) . Ax; 
err2=sum( (MTF_f it-MTF) .A2)  ; 
if  err2<min 
min=err2 ; 

ER= [ ER  err2 ] ; 
m= [x] ; 
end 

end 


%  plots  the  results 
f l=linspace (0,0.55,55)  ; 

MTF1= (sine (fl/ ( . 55) ) ) . Am(l) ; 

figured)  %  this  is  figure  4.4  in  chapter  IV 

plot ( f , MTF , 'o' , fl,MTFl, 'm' ) 

grid 

xlabel (' Frequency  in  cycles /miliradians ' ) 

ylabel ( 'MTF' ) 

axis ([0  0.6  0  1.2]) 


Naval  Postgraduate  School  -  CA  **************** 


%********* 
%  Type 
%  Name 


Thesis  Research 
Procedure 
mrtd__plot 

calculates  and  plots  MRTD 
Jul  /  1999 
1.0 

Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 
%  Routine  Information: 

%  1)  calculates  MRTD  (unpolarized  case) 

%  2)  calculates  MRTD  (polarized  case) 

%  3)  plots  the  results 


%  Date 
%  Version 
%  Author 
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%clear  all 
%  AGA  780  data 
SNRT=6 . 5 ; 
deltax=l . 1 ; 
deltay=l . 1 ; 
(mRad)  ; 

L=7; 

bar; 


(unpolarized  case) 

%  use  2.5  if  only  50%Pd  is  required 
%  in- scan  detector  subtense  (mRad) ; 
%  cross-scan  detector  subtense  in 

%  length  to  width  ratio  of  the  MRTD 


te=0.1;  % 

Fr=6 .25; 

Nos=l;  % 

NSS=1;  % 

FOVx=(7*pi/180) *1000; 
FOVy=(7*pi/180) *1000; 
F=0 .20; 

tau0=0.6;  % 

Nd=  1 ; 

Nsc=0 . 75 ;  % 

D=0 . 055; 

Deestar=4 . 12866el0 ; 
(calculated  form  NET) 
dndT=6 . 27e-5 ; 
(relative  to  300°K) 


eye  integration  time  (s) ; 

%  frame  rate  (Hz) ; 
overscan  ratio; 
serial  scan  ratio; 

%  in-scan  f ield-of -view  (mRad) ; 

%  cross-scan  f ield-of-view  (mRad) ; 

%  focal  length  (m) ; 
transmission  of  the  optics; 

%  number  of  detectors; 
scan  efficiency; 

%  aperture  (m) 

%  detectivity  (cm-HzA (1/2) -watt (-1) 
%  derivative  of  Planck's  Equation 


%  NET  calculation 

NET_numerator=20*F* (FOVx*FOVy*Fr*Nos*Nss ) A(0.5); 
NET_denominator=tauO* ( (pi*Nd*Nsc) . A (0 . 5 ) ) * (DA2 ) *deltax*delt 
ay*Deestar*dndT; 

NET=NET_numerator/NET_denominator ; 


%  spatial  frequency 
f0=linspace (0,0.55, 400)  ; 

%  noise  filter  factor  calculation 

rb=0.335;  % 

rox= (1+ (2*f 0*rb) .A2)  .A(-0.5)  ; 


%  fitting  curve  to  allow  better  interpolation 
MTFfit= (sine (f0/ .55) ) .A1.17; 


%  MRTD  calculation 
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MRTD_3  0  0=2  *  SNRT*NET*  (1 .  /MTFfit)  .  *  (  (rox)  .A0.5).*((((fO).A2)* 
deltax*deltay/L) . A0 . 5) * ( te*Fr*Nos*Nss) . A-0 .5; 

%  polarized  case  (T=262) 

taup=0.85;  %  transmission  of  the  polarizer 

dndT=4 . 00e-5 ;  %  derivative  of  Planck's  Equation 

(relative  to  262°K) 

%  NET  calculation 

NET_numerator=20*F* (FOVx*FOVy*Fr*Nos*Nss) A (0.5); 
NET_denominator=taup*tauO* ( (pi*Nd*Nsc) . A (0 .5)  )  * (DA2) *deltax 
*deltay*Deestar*dndT; 

NET=NET_ntunerator/NET_denominator; 

%  MRTD  calculation 

MRTD_pol=2*SNRT*NET* (1. /MTFfit) . * ( (rox) . A0 . 5 ).*(((( fO ). A2 ) * 
deltax*deltay/L) .A0.5) * (te*Fr*Nos*Nss) .A-0.5; 

%  plots  the  results 
figure (3) 

plot (f0,MRTD_300, fO,MRTD_pol, ' : ' ) 
xlabel (' spatial  frequency  in  mrad' ) 
ylabel ( ' MRTD ' ) 
axis ([0  0.55  0  20]) 

legend ( 'MRTD  unpolarized’  ,  'MRTD  polarized' , 2) 
grid 

%  zoom 

figure (4) 

plot (f0,MRTD_300, f0,MRTD_pol, ' : ' ) 
xlabel (' spatial  frequency  in  mrad') 
ylabel ( 'MRTD' ) 
axis ([0.3  0.55  0  12]) 

legend ('MRTD  unpolarized',  'MRTD  polarized ', 2 ) 
grid 


c^-k-k'k±'k*'kmk'k 
5 *  -k  *  *  *  *  *  * 

%  Type  : 

%  Name 
%  Function: 
%  Date  : 


Naval  Postgraduate  School  -  CA 
Thesis  Research 
Procedure 
sim_mrt2 

simulates  mrtd  measurements 
May  /  1999 
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%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  creates  a  four  bar  pattern 

%  2)  creates  a  simulated  output  for  regions  A,  B,  C,  ans 

%  D,  as  explained  in  figure  5.1  in  chapter  V. 

%********************************************************** 
%********************************************************** 

clear 

%  creates  a  four-bar  pattern 

m=zeros (250, 256) ; 

factor=8; 

a=factor*7 ; 

bar=ones (a, factor) ; 

m(  12 8- factor* 3 . 5  +  1 : 128+factor*3 . 5 , 12 5- factor *3 .5+1:125- 

factor*3 . 5+factor) =bar ; 

m( 12 8 -factor* 3 . 5+1 : 12 8+ factor* 3 .5,125- 

factor*3 . 5+l+2*factor : 125-factor*3 . 5+2*factor+factor) =bar; 
m( 12 8- factor *3 . 5+1 : 12 8+ factor* 3 .5,125- 

factor*3 . 5+l+4*factor : 125-factor*3 . 5+4*f actor+factor) =bar; 
m( 12 8 -factor* 3 . 5+1 : 128+factor*3 .5,125- 

factor*3 . 5+l+6*factor : 125-factor*3 . 5+6*factor+factor) =bar; 

%  calculates  the  pattern  fft 
f2=f ft2 (m) ; 
f2=f f tshif t ( f 2 ) ; 


a=[9  14  15  19];  %  scales  the  model  to  have  outputs  in 

regions  A,  B,  C,  and  D 
for  y=l : 4 

b=ones  (a  (y)  ,  a  (y)  )  ; 
for  i=l:a(y) 
for  j=l:a(y) 

b(i, j)=(-iA2-jA2) ; 
end 

end 
bl  =b  ; 

b2=f liplr (b) ; 
b3=f lipud (b) ; 
b4=f lipud (b2 ) ; 

BB=[b4  b3;b2  bl] ; 
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BB=BB+abs (BB(1)  )  ; 

BB=BB/inax  (max  (BB)  )  ; 

B=zeros(250,256); 

B(125-a(y)+l:125+a(y) , 128-a (y) +1 : 128+a (y) )=BB; 

%  multiplies  pattern  FFT  and  modeled  MTF 
test=f 2 . *B; 
test2=ifft2 (test)  ; 

%  plots  the  results 

figured)  %  this  is  figure  5.3  in  chapter  V 

subplot (2 , 2,y) 

surf (abs ( test2 ) ) 

axis( [100,150,100,150,  0  1]) 

axis  off 

end 


^•k‘k'k-k-k-k'k*'k~k**'k*'k-k'k*'k'k'k> 


% 


★  ★■A-****** 

Type 

Name 

Function 

Date 

Version 

Author 


% 

9- 


% 

% 


Naval  Postgraduate  School  -  CA  **************** 
Thesis  Research  ************** 

Function 
dl2temp 

converts  digital  level  to  temperature 
Jul  /  1999 
1.0' 

Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 
Routine  Information: 

1)  inputs: 

a)  digital  level 

b)  thermal  level 

c)  thermal  range 

d)  polarization  filter  number 

2 )  output : 

a)  temperature  in  degrees  centigrades 
,.********************************************************** 
U***'A,*^*‘A-*^*A-*A-*V:A:-A-ArA-*A-**^^A-Vr*A-*^*^***^*^****^****^':*:':?r'*'‘:*:'*‘:,r*'5,r 


function  tempC=dl2temp (dl , tl , tr , p) 
if  dl<=4095  &  dl>=0 
if  tl>=0  Sc  tl<=1000 

if  tr==2  |  tr==5  |  tr==10  |  tr==20  |  tr==50  |  tr==100 
|  tr==2  00  |  tr==500  |  tr==1000 
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%  transforms  digital  level  in  isothermal  units 
iu=dl* tr/4095+tl-tr/2 ; 

%  transforms  isothermal  units  in  temperature 
switch  p 
case  0 

A=2318 ;  B=1218;  C=0.742; 
case  4 

A=  5.679;  B=57.46;  C=0.918; 
case  5 

A=5 . 679 ;  B=57.46;  C=0.918; 

end 

tempK=B/ (log ( (A+iu) / (C*iu) ) ) ; 

%temperature  in  Kelvin 

tempC=tempK-273 . 15 ;  %temperature  in  Celcius 
else 

error (' Thermal  Range  must  be  2,  5,  10,  20,  50, 

100,  200,  500,  or  1000' ) 
end 
else 

error (' Thermal  Level  must  be  a  integer  between  0  and 

1000' ) 
end 
else 

error (' Digital  level  must  be  a  integer  between  0  and  4095') 
end 

if  tl<tr/2 

warning (' This  combination  of  thermal  range  and  thermal 
Level  may  result  in  complex  valued  calculations ' ) 
end 


%  *  *  *  k  *  -k  *  ★ 

%  Type  : 

%  Name  : 

%  Function: 
(filter  0) 


Naval  Postgraduate  School  -  CA  **************** 
Thesis  Research  ************** 

Procedure 
plotcal 

compares  calibrated  and  uncalibrated  data 


%  Date  :  Jul  /  1999 
%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt .  (BAF) 
%  Routine  Information: 

%  1)  plots  the  data  collected  from  the  AGA 
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% 


2)  fits  the  data 
%  3)  plots  the  calibrated  data 

%  4)  plots  the  calculated  MRTD 

£********************************************************** 
%* ******************************************* ************** 


clear 

%  measrued  data 

f = [ 0 . 1  0.2  0.3  0.3  0.4  0.4  0.43  0.43  0.43  0.45  0.47  0.48 
0.48  0.49  0.5  0.52  0.52] ; 

MRT= [0.1  0.3  1.9  0.7  2  1.3  0.4  10.9  10.5  1.5  2.1  35.4  17.9 
5.1  16.6  11.6  12.6] ; 

%  calibrated  data  (using  procedure  proposed  in  chapter  IV) 
MRTcal= [0.076277651  0.192184497  0.465116279  0.568643379 
1.958863859  1.352757544  0.539811066  3.367315416  3.37403599 
2.396166134  5.540897098  6.307911618  4.341498909  5.290456432 
6.566455696  7.16934487  9.01932713]; 

%  fitting  process 
%net=0 . 12 ; 

%das=l.l; 

%SL=0 . 075*net ; 

%SC=1 . 3  *net*das ; 

%ER=1; 

min=iel0; 

ER= [ ] ; 

for  B=0 . 6 : 0 . 01 : 0 . 8 

for  C=10: 0.01:12 

MRTD=0 . 05+B*f . *exp (C*f . A2 ) ; 
err2=sum( (MRTD-MRTcal) .A2) ; 
if  err2<min 
min=err2 ; 

ER= [ER  err2 ] ; 
m= [B  C] ; 
end 

end 

end 

min=lel0; 
for  Bl=0 : 0 . 1 : 4 

for  Cl=0 : 0.1:20 

MRTD1=0 . 05+Bl*f . *exp (Cl*f.A2) ; 
err2=STim(  (MRTDl-MRT)  .  A2)  ; 
if  err2<min 


min=err2 ; 

ER=[ER  err2 ] ; 
ml=[Bl  Cl] ; 

end 

end 

end 

%  spatial  frequency 
f 1=0: 0.01: 0.52; 

%  result  of  the  fitting  process 
MRTD=0 . 05+m(l) *fl . *exp (m (2 ) *f 1 . ^2 )  ; 

MRTD1=0 . 05+ml (1) *f 1 . *exp (ml (2 ) *f 1 . *2 ) ; 

%  gets  calculated  MRTD  from  mrtd_plot  function 
mrtd _plot ; 
close  all 

%  plots  the  results 

plot (fl, MRTD, f , MRTcal / 'x' , f 1 , MRTDl , ' : ' , f , MRT , 'o' , f0,MRTD_30 
0 , 

legend ( ' fitting  calibrated', 'calibrated', 'fitting 
uncalibrated' , 'uncalibrated' , 'calculated  MRTD' ,2) 
grid 

axis ( (0  0.55  0  40] ) 

xlabel (' spatial  frequency  in  cycles/mrad' ) 
ylabel ( ’ MRTD  -  °C ’ ) 


%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  plotcal 

%  Function:  compares  calibrated  and  uncalibrated  data 
%  (filter  5) 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  plots  the  data  collected  from  the  AGA 

%  2)  fits  the  data 

%  3 )  plots  the  calibrated  data 

o 
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,£********************************************************** 

clear 

%  measured  data 

f=[0.1  0.2  0.3  0.4  0.43  0.46  0.5]; 

MRT=[0.12  0.34  0.84  2.45  3.56  6.56  14.16]; 

MRTcal= [0. 07138608  0.214105793  0.717335611  2.437810945 
2.51590106  5.324675325  13.69439072]; 

%fitting  process 
%net=0 . 12 ; 

%das=l . 1; 

%SL=0 . 075*net; 

%SC=1 . 3  *net*das ; 

%ER=1 ; 
min=lel0; 

ER=  [  ]  ; 

for  B=0 .5:0.01:0.8 

for  C=10:0. 01:30 

%MRTD=SL+ ((SC*2*f)/2.646) . *exp (1.571* (ER^2 )*(f.A2)); 

MRTD=0 . 05+B*f . *exp (C*f . A2 ) ; 
err2=sum ( (MRTD-MRTcal ) .^2) ; 
if  err2<min 
min=err2 ; 

ER=[ERerr2]; 
m= [B  C] ; 
end 

end 

end 

min=lel0 ; 
for  Bl=0 : 0 . 1 : 4 

for  Cl— 0 : 0 . 1 : 40 

%MRTD=SL+ ( ( SC*2*f ) /2 . 646) . *exp (1 .571* (ER^2 )*(f.A2)); 
MRTD1=0.05+Bl*f  .*exp(Cl*f  ./'2)  ; 
err2=sum ( (MRTDl-MRT) .A2) ; 
if  err2<min 
min=err2 ; 

ER= [ER  err2] ; 
ml=[Bl  Cl] ; 
end 

end 

end 
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%  spatial  frequency 
f 1=0: 0.01: 0.52; 

%  result  of  the  fitting  process 
MRTD=0 . 05+m ( 1 ) *  f 1 . *exp (m  (2 )  * f  1 .  A2 ) ; 

MRTD1=0 . 05+ml (1) *f 1 . *exp (ml (2 ) *f 1 .  A2 )  ; 

%  plots  the  results 

plot (fl,MRTD, ' : ' , f , MRTcal , 'x' , f 1 , MRTD1 , f , MRT , 'o' ) 
grid 

axis ( [0  0.55  0  40] ) 

legend ( ' fitting  calibrated', 'calibrated', 'fitting 
uncalibrated' , 'uncalibrated' ,2) 
xlabel (' spatial  frequency  in  cycles /mrad ' ) 
ylabel ( 'MRTD  -  °C') 


%*********  Naval  Postgraduate  School  -  CA  **************** 
*  *  *  *  *  Thesis  Resecurch  ***★*****■*■★*** 

%  Type  :  Procedure 

%  Name  :  fig5_8a 

%  Function:  plots  calculated  vs  measured  MRTD  values  for 
%  filter  0  (unpolarized  case) 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes ,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  calculates  MRTD 

%  2)  plots  calculated  MRTD  vs  measured  MRTD 

^********************************************************** 
%* ********************************************************* 


%  measured  data  (unpolarized) 
clear 


%  first  set 

f_l= [0.3  0.4  0.43  0.48] ; 

MRT_1= [1.9  2  10.5  17.5] ; 

MRT_lcal= [0.51  1.65  3.20  3.91]; 

%  second  set 

f_2=[0.1  0.2  0.3  0.4  0.5]; 

MRT_2  =[0.06  0.28  0.67  1.28  16.56]; 
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MRT_2cal= [0.05  0.17  0.62  1.17  9.07]; 


%  third  set 

f_3= [0 . 1  0.2  0.3  0.4  0.5  0.52]; 

MRT_3  = [ 0 . 1  0.3  0.64  1.84  7.17  13.72]; 
MRT_3cal= [0.05  0.17  0.63  2.27  4.86  10.13]; 


%  AGA  780  data 
SNRT=  6.5; 
deltax=l . 1 ; 
deltay=l . 1; 

L=7; 
te=0 . 1; 

Fr=6 . 25 ; 

Nos=l; 

Nss=l ; 

FOVx= ( 7  *pi / 1 8  0 ) 
FOVy= (7*pi/180) 
F=0 .20; 
tau0=0.6; 

Nd=l ; 

Nsc=0 . 75 ; 

D=  0.055; 


(unpolarized  case) 

%  use  2.5  if  only  50%Pd  is  required  ; 

%  in-scan  detector  subtense  (mRad) ; 

%  cross-scan  detector  subtense  in  (mRad) ; 
%  length  to  width  ratio  of  the  MRTD  bar; 

%  eye  integration  time  (s); 

%  frame  rate  (Hz) ; 

%  overscan  ratio; 

%  serial  scan  ratio; 

*1000;  %  in-scan  f ield-of-view  (mRad) ; 

*1000;  %  cross-scan  f ield-of-view  (mRad) ; 

%  focal  length  (m) ; 

%  transmission  of  the  optics; 

%  number  of  detectors; 

%  scan  efficiency; 

%  aperture  (m) 


%  detectivity  (cm-HzA (1/2) -watt (-1)  =  (calculated  form  NET) 
Deestar=4 . 12866el0 ; 

%  derivative  of  Planck's  Equation  (relative  to  300°K) 
dndT=6 . 27e-5 ; 


%  NET  calculation 

NET_numer a t or = 2  0  *  F  * (FOVx*FOVy*Fr*Nos*Nss) A (0.5) ; 
NET_denominator= tauO *  (  (pi*Nd*Nsc)  .  A  ( 0 . 5) )  *  (DA2 )  *deltax*... 

deltay*Deestar*dndT; 

NET=NET_numerator/NET_denominator; 


%  spatial  frequency 
f 0=1 inspace (0, 0 . 55,400)  ; 

%  noise  filter  factor  calculation 

rb=0.335;  % 

rox= ( 1+ (2*f 0*rb) .A2) ,A(-0.5) ; 

%  fitting  curve  to  allow  better  interpolation 
MTFf it= (sine (f 0/ . 55) ) . A1 . 17 ; 
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MRTD_3  00=2*  SNRT  *  NET  *  (1.  /MTFfit)  .  *  (  (rox)  .  A0 . 5  ).*((((  fO )  .  A2  ) ... 
*deltax*deltay/L) . A0 . 5) * ( te*Fr*Nos*Nss) . A-0 . 5 ; 

%  plots  the  results 
figure (1) 

plot  ( f 1 ,  MRT_1 ,  '  x '  ,  f_2 ,  MRT_2  ,  's',  f_3  ,  MRT_3  ,  ... 

' o' , f 0 , MRTD_3 0 0 ,  ' -  '  ) 

legend  (’ First  set Second  set Third  set',  'calculated... 
MRTD ' ,2) 

grid 

axis ( [ 0  0.55  0  20] ) 

xlabel (' spatial  frequency  in  cycles/mrad' ) 
y label ( 'MRTD  -  °C' ) 


o 


Naval  Postgraduate  School 
Thesis  Research 


%  Type  :  Procedure 

%  Name  :  fig5_9 

%  Function:  plots  calculated  vs  measured  MRTD  values  for 
%  filter  4 


%  Date  :  Jul  /  1999 
%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 
%  Routine  Information: 


%  1)  calculates  MRTD 

%  2 )  plots  calculated  MRTD  vs  measured  MRTD 


%  measured  data  (polarized) 
clear 


f_l=[0.10  0.20  0.30  0.40  0.43  0.46  0.50  0.52]; 

MRT_1= [0.5600  0.7200  1.1700  2.0600  3.6700  6.4400  13.7200 
40.9500] ; 

MRT_lcal= [0.0974  0.1793  0.7638  1.1012  1.6141  2.6747  5.7369 
23.0554] ; 


%  AGA  780  data  (unpolarized  case) 

SNRT=6.5;  %  use  2.5  if  only  50%Pd  is  required  ; 

deltax=l.l;  %  in-scan  detector  subtense  (rnRad) ; 
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deltay=l . 1; 

L=7; 
te=0 . 1 ; 

Fr=6 .  25  ; 

Nos=l; 

Nss=l ; 

FOVx= (7*pi/180) *1000; 
FOVy= (7*pi/180) *1000; 


%  cross-scan  detector  subtense  in  (mRad) ; 
%  length  to  width  ratio  of  the  MRTD  bar; 

%  eye  integration  time  (s)  ; 

%  frame  rate  (Hz); 

%  overscan  ratio; 

%  serial  scan  ratio; 

%  in-scan  f ield-of-view  (mRad) ; 

%  cross-scan  f ield-of-view  (mRad) ; 


F=0.20;  %  focal  length  (m); 

tau0=0.6;  %  transmission  of  the  optics ; 

Nd=l;  %  niomber  of  detectors; 

Nsc=0 . 75 ;  %  scan  efficiency; 

D=  0.055;  %  aperture  (m) 

%  detectivity  (cm-HzA (1/2) -watt (-1)  =  (calculated  form  NET) 
Deestar=4 . 12866el0 ; 

taup=0 . 85 ;  %  transmission  of  the  polarizer 

%  derivative  of  Planck's  Equation  (relative  to  262°K) 
dndT=4 . 00e-5 ; 


%  NET  calculation 

NET_numer a t or =2  0  *  F  * (FOVx*FOVy*Fr*Nos*Nss ) A ( 0 . 5 ) ; 
NET_denominator=taup*tauO* ( (pi*Nd*Nsc) . A (0.5) ) * (DA2) *deltax 

*deltay*Deestar*dndT; 

NET=NET_numerator/NET_denominator; 

%  spatial  frequency 
f0=linspace (0, 0.55,400)  ; 

%  noise  filter  factor  calculation 

rb=  0.335;  % 

rox= (1+ (2*f 0*rb) . A2 ) . A (-0 . 5 )  ; 

%  fitting  curve  to  allow  better  interpolation 
MTFfit= (sine ( f 0 / . 55) ) .A1.17; 

%  MRTD  (polarized  case) 

MRTD_pol=2*SNRT*NET* (1. /MTFfit) . * ( (rox) . A0 . 5 ).*(((( f0) . A2 ) * 
deltax*deltay/L) . A0 . 5) * (te*Fr*Nos*Nss) .A-0.5; 

%  plots  the  results 
figure (1) 

plot  ( f _ 1 ,  MRT_1 ,  '  x '  ,  f  0  ,  MRTD__pol ,  '  -  '  ) 

legend (' Measured  values',  'calculated  MRTD  (polarized  - 
filter  4) ' , 2 ) 

grid 
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axis ( [0  0.55  0  50] ) 

xlabel (' spatial  frequency  in  cycles/mrad' ) 
ylabel ( * MRTD  -  °C 1 ) 

%********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  fig5_10 

%  Function:  plots  calculated  vs  measured  MRTD  values  for 
%  filter  5 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  calculates  MRTD 

%  2)  plots  calculated  MRTD  vs  measured  MRTD 

Sfei?****************************************************'*:**** 


%  measured  data  (polarized) 
clear 


f_l=[0.10  0.20  0.30  0.40  0.43  0.46  0.50  0.52]; 

MRT_1 =[0.1100  0.3400  0.8400  2.4500  3.5600  6.5600  14.1700 
22.3300] ; 

MRT_lcal= [0.0259  0.1578  0.4508  1.7640  1.5273  3.4273  5.7949 
9.4493] ; 


%  AGA  780  data 

SNRT=  6.5; 
deltax=l.l; 
deltay=l . 1 ; 

Ij=7  ; 
te=0 . 1; 

Fr=6 .25; 

Nos=l; 

Nss=l ; 

FOVx= (7*pi/180 ) 
FOVy= (7*pi/180 ) 
F=0 .20; 
tau0=0 . 6 ; 

Nd=  1 ; 

Nsc=0 . 75 ; 

D=0 .055; 


(polarized  case) 

%  use  2.5  if  only  50%Pd  is  required  ; 

%  in-scan  detector  subtense  (mRad) ; 

%  cross-scan  detector  subtense  in  (mRad) ; 
%  length  to  width  ratio  of  the  MRTD  bar; 

%  eye  integration  time  (s) ; 

%  frame  rate  (Hz) ; 

%  overscan  ratio; 

%  serial  scan  ratio; 

*1000;  %  in-scan  f ield-of-view  (mRad) ; 

*1000;  %  cross-scan  f ield-of-view  (mRad) ; 

%  focal  length  (m) ; 

%  transmission  of  the  optics; 

%  number  of  detectors; 

%  scan  efficiency; 

%  aperture  (m) 
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%  detectivity  (cm-HzA (1/2) -watt (-1)  =  (calculated  form  NET) 
Deestar=4 . 12866el0 ; 

taup=0.85;  %  transmission  of  the  polarizer 

%  derivative  of  Planck’s  Equation  (relative  to  262°K) 
dndT=4 . 00e-5 ; 

%  NET  calculation 

NET_numerator=2 0  *F*  (FOVx*FOVy*Fr*Nos*Nss)  A  ( 0 . 5 )  ; 
NET_denominator=taup*tauO* (  (pi*Nd*Nsc)  .A(0.5)  )  *  (DA2 )  *deltax 

*deltay*Deestar*dndT; 

NET=NET_numerator/NET_denominator; 

%  spatial  frequency 
f0=linspace(0,0.55,400)  ; 

%  noise  filter  factor  calculation 

rb=0 . 335 ;  % 

rox= (l+(2*f0*rb)  .A2) .A(-0.5)  ; 

%  fitting  curve  to  allow  better  interpolation 

MTFf it= (sine ( f 0/ . 55 ) ) . A1 . 17 ; 

%  MRTD  (polarized  case) 

MRTD_pol=2*SNRT*NET* (1. /MTFfit)  . * ( (rox)  . A0 .5)  .* ( ( (  (f0)  . A2) * 
deltax*deltay/L) . A0 . 5) * ( te*Fr*Nos*Nss) .A-0.5; 

%  Plots  the  results 
figured) 

plo  t  ( f _ 1 ,  MRT_l,'x',  f  0  ,  MRTD__pol ,  '  -  '  ) 

legend (' Measured  values',  'calculated  MRTD  (polarized  - 
filter  5 ) ' ,2) 

grid 

axis ( [0  0.55  0  50] ) 

xlabel (' spatial  frequency  in  cycles/mrad' ) 
ylabel ( 'MRTD  -  °C' ) 


-k  *  *  *  ~k  *  *  * 

%  Type  : 

%  Name  : 

%  Function: 
%  Date  : 


Naval  Postgraduate  School 
Thesis  Research 
Procedure 


dndt 

plots  Planck's  equation  and  its  derivative 
Jul  /  1999 
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%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  calculates  M  (Planck's  equation) 

%  2)  calculates  N  (radiance) 

%  3)  calculates  the  derivative  of  Planck's  equation 

%  4)  calculates  the  derivative  of  N 

%  5)  Plots  M,  dMdt,  and  dNdt 

^********************************************************** 


clear 
close  all 
lambda=8 : 0 . 001 : 12 ; 

%  inputs  the  temperature 
T=input (' Enter  T :>>'); 

%  Planck's  Equation  Constants 
cl=3 .7418e4; 
c2=l . 4388e4; 

%  Planck's  Equation 
Ml= (cl . / (lambda. *5) )  ; 

M2=exp (c2 . / (lambda*T) ) -1; 

M=Ml . /M2 ; 

Ni=mean (M) /pi 

%  The  derivative  of  Planck's  Equation 

dMdT=M.  *c2  .  *exp  (c2  .  /  (lambda*T)  )  .  /  (lambda*  (T/'2)  .  *exp  (c2  .  /  (la 

mbda*T) ) ) ; 

dNdT= dMdT /pi; 

dNdT=0 . 001*cumsum(dNdT)  ; 

%  Plot  the  results 
figure (1) 
subplot (2,1,1) 
plot ( lambda, M) 
grid  on 

xlabel ( ' lambda  -  micrometer' ) 

ylabel ( 'M' ) 

hold  on 

subplot (2,1,2) 

plot ( lambda, dMdT) 

grid  on 
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xlabel ( ' lambda  -  micrometer' ) 
ylabel ( 'dMdT' ) 
figure  (2 ) 

plot (lambda,  dNdT) 
grid 

dndt=dNdT( length (dNdT) ) 
orient  tall 


9.*  ********************************************************* 

**********  Naval  postgraduate  School  -  CA  **************** 

Thesis  Research  ************** 

Procedure 
enne 

calculates  radiance  for  a  set  of  temperatures 
Jul  /  1999 
1.0 

Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  calculates  radiance 

%  2)  plots  the  result 

%  *  the  result  is  presented  in  Appendix  D 

^* **********************  *********************************** 

%********************************************************** 


Qr-k'k'k-k'k'k'k** 

%  Type 
%  Name 
%  Function 
%  Date 
%  Version 
%  Author 


format  long 
clear 
close  all 
lambda=8 : 0 . 001 : 12  ; 

Ni=  [  ]  ; 

%  initializes  the  sset  of  temperatures 
for  T=240 : 320 ; 
cl=3 .7418e4; 
c2=l . 4388e4 ; 

%  calculates  Planck's  Equation 
Ml= (cl. / (lambda. A5) ) ; 

M2=exp(c2 . / ( lambda* T) ) -1; 

M=M1 . /M2 ; 

Ni= [Ni  mean (M) /pi]; 

%  outputs  string  of  radiance  values 
fprintf ( 1 , [ ' %6 . Of  %10.8e\n'] , T, mean (M) /pi ) ; 
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end 


%  plots  the  results 
plot (240 : 320 , Ni) 
grid 

xlabel ( 'T' ) 
ylabel ( 'N' ) 

axis ( [240  320  2e-4  1.4e-3]) 


*********  Naval  Postgraduate  School  -  CA  **************** 


% 

% 


'6 

O. 


% 


Thesis  Research 
Type  :  Procedure 

Name  :  rangeimp4 

Function:  estimates  range  improvement  using  polarization 
degree 

Date  :  Jul  /  1999 

Version  :  1.0 

Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt .  (BAF) 
Routine  Information: 

1)  calculates  the  MRTD  as  a  function  of  range 

2)  plots  target- to-background  temperature  difference 
for  four  polarization  cases  (0%,  10%,  20%,  and  30%); 
using  table  in  Appendix  D 


format  long 
clear 
close  all 
lambda=8 : 0 . 001 : 12 ; 

%  results  presented  in  table  5.2  (using  table  D.l) 
deltaTl=12 . 7 ; 
deltaT2=18; 
deltaT3=23 .7; 

%  Using  Searad  data  with  the  following  parameters: 

%  *  Lowtran7  was  selected 

%  •  Type  of  atmospheric  path:  Slant  path  between  two 
%  altitudes 

%  •  Surface  albedo  of  earth:  assumed  blackbody 
%  •  Navy  maritime  model 
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%  •  No  clouds  or  rain 
%  •  Altitude  1:  0 
%  •  Altitude  2:  300ft 
%  •  Wind  speed:  7.5  m/s 
%  •  Wind  speed:  7.5  m/s  (24h  average) 

%  •  Relative  Humidity: 

%  •  Airmass  character:  3 

%  •  Visibility:  entered  0  (means  that  the  program  will  use 
%  the  standard  data  for  the  chosen  condition 

%  <Mid- latitude  summer >) 

Range=linspace(0, 40,400) ; 

%  Searad  Data 
Title= [ ' 300ft 1 ] ; 

r=[0  1  2  3  4  5  10  15  20  25  30  35  40]; 

t=[l  0.7326  0.5784  0.4634  0.3745  0.3044  0.1155  0.0465 
0.0194  0.0083  0.0036  0.0016  0] ; 
y=  [  ]  ; 

%  linearly  filling  the  spaces  between  two  data  points 
for  j=l:12 

temp= 1 inspace (t (j ) , t ( j+1) , (r(j+l)-r(j) ) *10+1) ; 
temp ( length ( temp) ) = [ ] ; 
y=  [y  temp] ; 
end 

%  This  function  calculates  the  MRTD 

mrtd_plot 

close  all 

R0=f 0*2*15;  %  15m  =  critical  dimension  of  ship 

%  Plotting 
figure (1) 

plot (Range, y* 10, Range,  y*deltaTl, Range,  y*deltaT2 , Range, 

y*deltaT3 , R0 , MRTD_3  00 , R0 , MRTD_pol ) 

xlabel ( 'Range  (km) 1 ) 

ylabel ( ' TBTD ' ) 

axis ( [0  25  0  5]  ) 

grid 

figure (2) 
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plot (Range, y*10 , Range,  y*deltaTl, Range,  y*deltaT2 , Range, 

y*deltaT3 ,R0,MRTD_300,R0,MRTD __pol) 

xlabel ( 'Range  (km) ' ) 

ylabel ( ' TBTD ' ) 

axis ( [8  14  0.5  2.5] ) 

grid 

gtext ( 'POL=0%' ) 

gtext ( ' POL=10% ' ) 

gtext ( 'POL=20%' ) 

gtext ( 'POL=30%' ) 

gtext ('MRTD  -  with  polarizer') 

gtext ('MRTD  -  without  polarizer') 

figure (3 ) 

plot ( R0 , MRTD_3  0 0 , R0 , MRTD_pol ,  ' :  '  ) 

axis ( [0  16  0  10] ) 

xlabel (' Range  (km)') 

ylabel ( 'MRTD' ) 

grid 

legend ('MRTD  -  without  polarizer' , 'MRTD  -  with 
polarizer ' , 0 ) 


^********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  range imp 5 

%. Function:  estimates  range  improvement  using  Lagaras 
%  [Ref.  11]  data 

%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt .  (BAF) 

%  Routine  Information: 

%  1)  calculates  the  MRTD  as  a  function  of  range 

%  2)  plots  target-to-background  temperature  difference 

%  for  polarized  and  unpolarized  cases 

%********************************************************** 
%********************************************************** 


format  long 
clear 
close  all 
lambda=8 : 0.001:12; 
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%  Using  Searad  data  with  the  following  parameters: 

%  •  Lowtran7  was  selected 

%  •  Type  of  atmospheric  path:  Slant  path  between  two 
%  altitudes 

%  •  Surface  albedo  of  earth:  assumed  blackbody 
%  •  Navy  maritime  model 
%  •  No  clouds  or  rain 
%  •  Altitude  1:  0 

%  •  Altitude  2:  300ft 

%  •  Wind  speed:  7.5  m/s 

%  •  Wind  speed:  7.5  m/s  (24h  average) 

%  •  Relative  Humidity: 

%  •  Airmass  character:  3 

%  •  Visibility:  entered  0  (means  that  the  program  will  use 
%  the  standard  data  for  the  chosen  conditions 

%  <Mid-latitude  summer>) 

Range=l inspace (0,40, 400) ; 

%  Lagaras  data 

rangel= [1 . 5  3.5  5  7  8  9  11  13  15.5  16.5  19  20  21  22  25 

25.5  27  29  31  33  40] ; 

unp= [15.7  12.8  11  9  8  7.5  6.4  5  4  3.5  3  2.6  2.5  2.4  2  1.8 

1.5  1.4  1  0.7  0 ] ; 

range 2= [1.5  3.5  4.5  7  8  9  11  13  15 . 5  16 . 5  19  20  22  25  40] ; 

pol= [36  29  27  21  19  17.5  14.5  12  8.9  8.15  6  5.5  4.6  3.5 

0]; 

%  interpolation  process 
XI=0 : 0 . 1 : 40 ; 

YIl  =  INTERPl (r angel , unp, XI , 1  cubic  1 ) ; 

YI2  =  INTERPl (range2,pol, XI, ' cubic ') ; 

%  This  function  calculates  the  MRTD 

mrtd_plot 

close  all 

R0=f 0*2*15;  %  15m  =  critical  dimension  of  ship 

%  Plotting 

figure (1)  %  this  is  figure  5.15  in  chapter  V 

plot (rangel, unp,  'bo',  range2,pol,  ' rx' ,XI, YIl, 'b' , 

XI , YI2 , 'r ' ) 
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grid 

xlabel ( ' Range  (km) ' ) 
ylabel ( ' TBTD ' ) 

legend ( ' Unpolarized' , ' Polarized' ,  0 ) 

figure (2)  %  this  is  figure  5.16  in  chapter  V 

plot (XI, YI1, 'b' ,XI, YI2, 'r: ' , RO , MRTD_300 , 'b' , RO ,MRTD_pol , 'r: 

'  ) 

xlabel ( 'Range  (km) ' ) 
ylabel ( 'TBTD' ) 
axis ( [0  25  0  25] ) 
grid 

legend ( 'Unpolarized' , ' Polarized'  ,  0) 

figure (3)  %  this  is  figure  5.16  in  chapter  V 

plot  (XI,  YI1,  'b'  ,XI,  YI2,  'r:  '  ,  RO , MRTD_3 0 0  ,  'b'  ,  RO , MRTD__pol ,  'r: 

'  ) 

xlabel ( 'Range  (km) ' ) 
ylabel ( 'TBTD' ) 
axis ( [13  16  3  10] ) 
grid 

legend ( 'Unpolarized' , 'Polarized' , 0) 


^********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  fig6_ld 

%  Function:  compares  experimental  and  calculated  data 
%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  1)  plots  calculated  MRTD  vs  measured  MRTD 

%  (unpolarized) 

%  2)  plots  calculated  MRTD  vs  measured  MRTD  (polarized) 

%********************************************************** 
^********************************************************** 


clear 

%  measured  data  (unpolarized) 
f_0=[0.1  0.2  0.3  0.4  0.5  0.52]; 
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MRT_0= [0.1  0.3  0.64  1.84  7.17  13.72]; 

MRT_0cal= [0.05  0.17  0.63  2.27  4.86  10.13]; 

%  measured  data  (polarized  -  filter  4) 

f_4= [0.10  0.20  0.30  0.40  0.43  0.46  0.50  0.52]; 

MRT_4= [0 . 5600  0.7200  1.1700  2.0600  3.6700  6.4400  13.7200 
40.9500] ; 

MRT_4cal=[ 0.0974  0.1793  0.7638  1.1012  1.6141  2.6747  5.7369 
23.0554] ; 


%  measured  data  (polarized  -  filter  5) 

f_5= [0.10  0.20  0.30  0.40  0.43  0.46  0.500.52]; 

MRT_5= [0 . 1100  0.3400  0.8400  2.4500  3.5600  6.5600  14.1700 
22.3300]; 

MRT_5cal= [0 . 0259  0.1578  0.4508  1.7640  1.5273  3.4273  5.7949 
9.4493] ; 


%  AGA  780  data 

SNRT=6 . 5 ; 
deltax=l . 1; 
deltay=l . 1; 

L=7  ; 
te=0 . 1 ; 


(polarized  case) 

%  use  2.5  if  only  50%Pd  is  required  ; 

%  in-scan  detector  subtense  (mRad) ; 

%  cross-scan  detector  subtense  in  (mRad) ; 
%  length  to  width  ratio  of  the  MRTD  bar; 

%  eye  integration  time  (s) ; 


Fr=6 . 25 ;  %  frame  rate  (Hz)  ; 

Nos=l;  %  overscan  ratio; 

Nss=l ;  %  serial  scan  ratio; 

FOVx= (7*pi/180) *1000;  %  in-scan  f ield-of-view  (mRad); 

FOVy= (7*pi/180) *1000;  %  cross-scan  f ield-of-view  (mRad) ; 


F=0 .20; 
tau0=0 . 6 ; 
Nd=l  ; 
Nsc=0 . 75 ; 
D=  0 . 0  5  5 ; 


focal  length  (m) ; 
transmission  of  the  optics ; 
nimber  of  detectors; 
scan  efficiency; 
aperture  (m) 


%  detectivity  (cm-Hz^ ( 1/2 ) -watt ( -1)  =  (calculated  form  NET) 

Deestar=4 . 12866el0; 

taup=0 . 85 ;  %  transmission  of  the  polarizer 


%  derivative  of  Planck's  Equation  (relative- to  262°K) 
dndT=4 . 00e-5 ; 


%  NET  calculation 

NET_numerator=20*F* (FOVx*FOVy*Fr*Nos*Nss) ^ (0 . 5) ; 
NET_denominator=taup*tauO* ( (pi*Nd*Nsc) .^(0.5) ) * (D^2) *deltax 
*deltay*Deestar*dndT; 
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NET=NET_numerator/NET_denominator ; 

%  spatial  frequency 
f 0=linspace (0,0.55, 400)  ; 

%  noise  filter  factor  calculation 
rb=0 . 335 ;  % 

rox=(l+(2*f0*rb) .A2) .A(-0.5)  ; 

%  fitting  curve  to  allow  better  interpolation 
MTFf it= (sine (f0/ . 55) ) . A1 . 17; 

%  MRTD  (polarized  case) 

MRTD_pol=2*SNRT*NET* (l./MTFfit) .*((rox) .A0.5) .*((((f0) .A2)* 
deltax*deltay/L) . A0 . 5 ) * (te*Fr*Nos*Nss) .  A-0 . 5 ; 

%  AGA  780  data  (unpolarized  case) 

dndT=6 . 27e-5 ;  %  derivative  of  Planck's  Equation 

(relative  to  300°K) 

%  NET  calculation 

NET_numerator=2  0  *F* (FOVx*FOVy*Fr*Nos*Nss) A ( 0 . 5 ) ; 
NET_denominator=tauO* ( (pi*Nd*Nsc) . A ( 0 . 5 ) ) * (DA2 ) *deltax*delt 
ay*Deestar*dndT; 

NET=NET_numerator/NET_denoiuinator; 

%  MRTD  (unpolarized  case) 

MRTD_3  00=2*SNRT*NET* (1. /MTFf it)  . * ( (rox)  . A0 . 5 ).*(((  (fO )  . A2 ) * 
deltax*deltay/L) . A0 . 5) * (te*Fr*Nos*Nss) . A-0 . 5; 

%  plots  the  results 
figure (2 ) 

plot ( f_0 ,  MRT_0 , 1 bx 1 ,  f_4,  MRT_4 , ’ ro 1 , f _5 ,  MRT_5,'rs', 
f  0  ,  MRTD_3  00,'  b-  '  ,  f  0  ,  MRTD__pol ,  '  r  :  '  ) 

legend (' Measured  unpolarized ' , 'Measured  vertical  pol . 

' , 'Measured  horizontal  pol. ' ,  'calculated  MRTD  unpol . ' , 

'calculated  MRTD  pol. ',2) 

grid 

axis ( [0  0.55  0  20] ) 

xlabel (' spatial  frequency  in  cycles/mrad 1 ) 
ylabel ( 'MRTD  -  °C') 
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oV  oV 


********************************************************** 
*********  Naval  postgraduate  School  -  CA  **************** 
5*********  Thesis  Research  ************** 

5  Type  :  Procedure 

Name  :  readptw2 


% 

% 


Function:  reads  *.ptw  files  generated  by  PTRwin 


Date  :  May  /  1999 
%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 

%  Routine  Information: 

%  I)  It  reads  a  ’'.ptw  file  and  plots  it  as  a  figure  and 
%  also  plots  a  3D  view  of  the  image. 

%  2)  It  allows  the  user  to  select  a  detail  to  be 

%  analyzed. 

%  3)  It  plots  the  profile  of  the  detail. 

%  4)  It  allows  the  user  to  select  the  peaks  and  valleys 

%  of  the  profile. 

%  5)  It  performs  a  fourier  analysis  of  the  profile. 

%  6)  It  generates  a  string  output  to  the  file 

%  " summary . dat "  with  the  basic  information  of  the 

%  file:  name,  spatial  frequency,  target  temperature, 

%  background  temperature,  thermal  range,  thermal 

%  level,  peak  temperature,  valley  temperature,  type  of 

ct  filtsir  usscl . 

°******************'****************************************- 

„********************************************************** 


%  creates  output  file  'summary.dat' 
f id0=f open ( ' summary . dat ' , 'wr' ) ; 
f=l; 

while  f~=0 
%open  the  *.ptw  file 
f=uigetf ile ( ' * .ptw' ) ; 
fid=fopen (f , ' r ' )  ; 

pol=input ( 'Enter  the  filter  used  (0,  4,  5):'); 

%  locate  the  thermal  level  in  the  file 
gapl=6553  6 ; 

status=f seek ( fid, gapl , 'bof ' ) ; 
level=f read (fid, 1 , ' intl6 ' ) ; 

%  locate  the  thermal  range  in  the  file 
gap2=6553  6+40 ; 

status=f seek ( fid, gap2 , 'bof' ) ; 
range=f read (fid, 1 , ' intl6 ' ) ; 
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%  locate  the  image  in  the  file 

gap=65536+3072 ; 

status=f seek ( fid, gap, 'bof ' ) ; 

A=fread(fid, [16000, 4] , 'intl6' ) ; 
f ieldl=reshape (A( : , 1) , 250, 64) ; 
f ield2=reshape ( A ( : ,2) ,250, 64) ; 
f ield3=reshape (A( : , 3 ) , 250 , 64 ) ; 
f ield4=reshape (A( : , 4) ,250, 64) ; 
picture= [ ] ; 

%  build  the  image  with  the  fields 
for  n=l:64 

temp=  [f  ieldl  (  :  ,  n)  '  ;  f  ield2  ( :  ,n)  ' ,-  f  ield3  ( :  ,  n)  '  ;  f  ield4  ( :  ,  n)  '  ] 
picture= [picture; temp] ; 

end 

close  all 


%  plot  the  image  as  a  surface 

figured) 

surf (picture) 

view (0,90) 

%  select  a  detail 

title ('click  on  the  upper  left  corner  of  the  detail') 

[ dl , d2 ] =ginput ( 1 ) ; 

title ('click  on  the  lower  right  corner  of  the  detail') 
[d3 , d4 ] =g input ( 1 ) ; 
title ( ' ' ) 

%  plot  the  datail 
figure (2) 

detail=picture (round(d4) :round(d2) ,round(dl) :round(d3) ) ; 

%%%integers 

surf 1 (detail) 

shading  interp; 

colormap (pink) 

title (['3D  view  of  the  output  -  file:  ',f]) 

%  calculate  and  plot  the  profile 
figured) 

[x,y] =size (detail) ; 

prof ile=detail (round (x/2 ) , :) ;  %%%integers 
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%  transform  digital  level  in  temperature 
for  z=l : length (profile) 

profile (z) =dl2 temp (profile (z) , level, range, pol) ; 
end 

plot (prof ile) 
grid 

%  select  peaks  and  valleys 

title (' select  the  four  apparent  peaks,  using  the  mouse') 

[px, py] =ginput (4 ) ; 

title (' select  the  three  apparent  valleys,  using  the  mouse') 
[vx, vy] =ginput (3 ) ; 
peak=mean (py) ; 
valley=mean (vy) ; 

text (mean (px) , peak, [ 'mean  peak  is:  ' , num2str (peak) , 1 °C ' ] ) ; 
text (mean (px) , valley, [ 'mean  valley  is: 

' , num2str (valley) , ' °C ' ] ) ; 

title ( [ 'Horizontal  profile  of  the  output  -  file:  ’,f]) 
orient  tall 

%  allow  user  to  enter  file  parameters 
clc 

spatial_freq= input (' Enter  the  spatial  frequency  used  in 
this  file: ' ) ; 

Tb= input ( ' Enter  the  recorded  background  temperature 
(Tb) :'); 

Tt=input (' Enter  the  recorded  target  temperature  (Tt) : ' )  ; 
MRTD=Tb-Tt ; 

MRTDcal=0 . 1*MRTD/ (peak-valley) ; 

fprintf (fidO , [f,  '%6.2f  %6.2f  %6.2f  %6.2f  %6.2f  %8.4f  %8.4f 

%8.4f  %8.4f  %li\n ' ] , . . . 

spatial_f req, Tt , Tb, level , range , peak, valley , MRTD, MRTDcal , pol 

)  ; 

end 

f close ( ' all ' )  ; 
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'*■*•*•*■*■*•*** 
■  'k'k'k'k'k'krk'k'k 


************** 


Type 

Name 


% 


% 

o. 


% 

9: 


Naval  Postgraduate  School  -  CA 
Thesis  Research 
Procedure 
readmtf 

Function:  reads  MTF  information  from  selected  *.ptw  files 
generated  by  PTRwin. 

May  /  1999 
1.0 

Edson  Fernando  da  Costa  Guimaraes,  Capt .  (BAF) 


Date 
Version 
Author 

Routine  Information: 

1)  It  reads  a  *.ptw  file  and  plots  it  as  a  figure  and 
also  plots  a  3D  view  of  the  image. 

2)  It  allows  the  user  to  select  a  detail  to  be 
analyzed. 

3)  It  plots  the  profile  of  the  detail. 

4)  It  allows  the  user  to  select  the  peaks  and  valleys 
of  the  profile. 

5)  It  generates  a  string  output  with  the  basic 

%  information  of  the  file:  name,  peak  temperature,  and 

%  valley  temperature. 

96*  ********************************************************* 
Q.********************************************************** 


while  f~=0 

%  open  the  *.ptw  file 
clear 

f=uigetfile ( 'mtf * .ptw' ) ; 
f id=f open ( f , ' r ' ) ; 

%  locate  the  thermal  level  in  the  file 
gapl=65536; 

status=fseek(fid,gapl, 'bof ' ) ; 
level=f read (fid, 1 , ' intl6 ' ) ; 

%  locate  the  thermal  range  in  the  file 
gap2=65536+40 ; 

status=fseek(fid,gap2, 'bof' ) ; 
range=f read (fid, 1 , ' intl6 ' ) ; 

%  locate  the  image  in  the  file 
gap=65536+3072 ; 
status=f seek ( fid, gap, 'bof') ; 
A=fread(fid, [16000,4] , 'intl6' ) ; 


f ieldl=reshape (A( : , 1) , 250, 64) ; 
f ield2=reshape (A ( : , 2) , 250 , 64) ; 
field3=reshape ( A { : ,3) ,250, 64)  ; 
f ield4=reshape (A( : , 4) , 250 , 64) ; 
picture= [ ] ; 

%  build  the  image  with  the  fields 
for  n=l : 64 


temp= [fieldl ( : ,n) ' ; field2 ( : ,n) ' ; f ield3 ( : , n) ' ; f ield4 ( : , n) '] 
picture= [picture; temp] ; 
end 

close  all 

%  plot  the  image 
figure (1) 
imagesc (picture) ; 
title  (['Image:  ',f]); 

%  plot  the  image  as  a  surface 
figure (2) 
surf (picture) 
view (0,90) 

%  select  a  detail 

title ('click  on  the  upper  left  corner  of  the  detail') 

[ dl , d2 ] =ginput ( 1 ) ; 

title ('click  on  the  lower  right  corner  of  the  detail') 

[ d3 , d4 ] =ginput ( 1 ) ; 
title ( ' ' ) 

detail=picture (round (d4) :round(d2) ,round(dl) :round(d3) ) ; 

%%%integers 

%  calculate  and  plot  the  profile 
f igur e ( 4 ) 

[x,y] =size (detail) ; 

prof ile=detail ( round (x/2) , : ) ;  %%%integers 

%  select  peaks  and  valleys 

plot (profile) 

grid 

title (' select  the  four  apparent  peaks,  using  the  mouse') 
[px, py] =ginput ( 4 )  ; 
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title (' select  the  three  apparent  valleys,  using  the 
mouse' ) 

[vx, vy] =ginput (3 ) ; 
peak=mean (py) ; 
valley=mean (vy) ; 

text (mean (px) ,peak, ['mean  peak  is:  ' , num2str (peak) , ' 
DL']); 

text (mean (px) , valley, [ 'mean  valley  is: 

' , num2str (valley) , '  DL ' ] ) ; 

title ([ 'Horizontal  profile  of  the  output  -  file:  ',f]) 
orient  tall 

%  outputs  string  of  information 

fprintf (1, [f , '%8 . 4f  %8.4f\n' ] , peak, valley) ; 

end 


%********************************************************** 
%*********  Naval  Postgraduate  School  -  CA  **************** 
%*********  Thesis  Research  ************** 

%  Type  :  Procedure 

%  Name  :  filt 

%  Function:  set  filter  to  program  readgui 
%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes,  Capt.  (BAF) 
^* ********************************************************* 
<£********************************************************** 


function  fig  =  filt() 

%  This  is  the  machine-generated  representation  of  a  Handle 
%  Graphics  object  and  its  children.  Note  that  handle  values 
%  may  change  when  these  objects  are  re-created.  This  may 
%  cause  problems  with  any  callbacks  written  to  depend  on 
%  the  value  of  the  handle  at  the  time  the  object  was  saved. 
%  To  reopen  this  object,  just  type  the  name  of  the  M-file 
%  at  the  MATLAB  prompt.  The  M-file  and  its  associated  MAT 
%  -file  must  be  on  your  path. 

load  filt 

hO  =  figure ( 'Color ', [0 . 8  0.8  0.8],  ... 

' Colormap' ,mat0 ,  ... 
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' PointerShapeCData' ,matl,  . . . 

' Position' , [601  361  363  85],  ... 

'Tag' , 'Figl' ) ; 

hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units 'points ' ,  ... 

'Callback' , 'warning  off, global  p,p=0;warning  on, close', 
' ListboxTop' ,0 ,  ... 

'Position',  [23.58620689655173  6.827586206896553 
47.17241379310346  17.37931034482759],  ... 

' String' , ' 0 ' ,  ... 

' Tag ' , ' Pushbuttonl ' )  ; 
hi  =  uicontrol (' Parent ' ,h0,  ... 

'Units 'points ' ,  ... 

'Callback' , 'warning  off, global  p, p=5 ;warning  on, close', 
' ListboxTop' , 0 ,  ... 

'Position' , [162  6.827586206896553  47.17241379310346 
17.37931034482759],  ... 

' String' , ' 5 ' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 
hi  =  uicontrol ( ' Parent ' , hO ,  ... 

'Units ', 'points ' ,  ... 

'Callback' , 'warning  off, global  p, p=4 ;warning  on, close', 
'ListboxTop' , 0,  ... 

'Position', [94.34482758620692  6.206896551724139 
47.17241379310346  17.37931034482759],  ... 

' String' , ' 4 ' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 
hi  =  uicontrol ( ' Parent ' , hO ,  ... 

'Units' , 'points ' ,  ... 

'BackgroundColor ' , [0.752941176470588  0.752941176470588 
0.752941176470588],  ... 

'ListboxTop' , 0,  ... 

'Position',  [44.6896551724138  29.79310344827587 
141.5172413793104  15.51724137931035],  ... 

' String' ,' Select  Filter',  ... 

' Style ' , ' text ' ,  ... 

' Tag ' , ' StaticTextl ' ) ; 
if  nargout  >0,  fig  =  hO;  end 
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*********  Naval  Postgraduate  School  -  CA: **************** 


********* 


************** 


Thesis  Research 
Function 
guil 

support  the  Graphical  User  Interface  (GUI)  for 
program  readgui . 

Jul  /  1999 

1.0 

Edson  Fernando  da  Costa  Guimaraes,  Capt .  (BAF) 
Routine  Information: 

1)  loads  figure  from  file 

2)  gets  a  detail 

3)  calculates  the  profile 

%  4)  allows  selection  of  peaks  and  valleys 

%  5)  allows  frequency  domain  analysis 

%  6)  plots  detail  in  3D 

CU  ********************************************************** 


Type 

Name 

Function 

Date 

Version 

Author 


% 

% 


function  guil (action) 
global  p;  %  polarization 
colormap ( ' default ' ) 

global  fid  picture  f  detail  profile  level  range  vy  peak 
valley 

switch  action 
case  'load' 

f=uigetf ile ( ' * .ptw' )  ; 
f id=fopen ( f , ' r ' )  ; 
gapl=65536; 

s tatus=f seek (fid, gapl, 'bof ' ) 
level=fread(fid, 1, 'intl6' ) ; 
gap2=65536+40 ; 

status=f seek ( fid, gap2 , ' bof ' ) 
range=fread ( fid, 1 , ' intl6 ' ) ; 


%  loads  figure  from  file 
%  opens  file 

%  skips  the  header 

%  reads  thermal  level 


%  reads  thermal  range 


gap=65536+3072; 
status=f seek ( fid, gap, 'bof' ) ; 
A=fread(fid, [16000, 4] , ' intl6 ' ) ; 
f ieldl=reshape (A ( : , 1) ,250,64) ; 
f ield2=reshape (A ( : ,2) ,250, 64)  ; 
f ield3=reshape (A ( : , 3 ) ,250,64) ; 
f ield4=reshape (A( : , 4) , 250, 64) ; 
picture= [ ] ; 


%  reads  frame 
%  extracts  field  #1 
%  extracts  field  #2 
%  extracts  field  #3 
%  extracts  field  #4 
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%builds  the  picture  from  fields 
for  n=l : 64 

temp= [f ieldl  ( : ,n)  ' ;field2 ( : ,n) ' ;field3 ( : ,n) ' ; f ield4 ( : ,n) ' 

picture= [picture; temp] ; 
end 

%  plots  image 
imagesc (picture) ; 
title  (['Image:  ',f]); 


case  'detail'  •%  gets  a  detail 
imagesc (picture) ; 
title  (['Image:  ',f]); 

%  allows  selection  of  detail  limits 

title ('click  on  the  upper  left  corner  of  the  detail') 
[xl,x2]=ginput (1)  ; 

title ('click  on  the  lower  right  corner  of  the  detail') 

[x3 , x4 ] =ginput ( 1 ) ; 
title ( ' ' ) 

%  plots  detail 

detail=picture (min (round (x4) ,round(x2)) : max (round (x4) , rou 
nd  (x2 )  )  , min ( round (xl ) , round (x3 ) ) : max ( round (xl ) , round (x3 ) ) ) ; 
%%% integers 

surf (detail) 

[a,b]=size(detail)  ; 

axis  (  [1  b  1  a  0  max  (max  (detail)  )  ]  ) 
view (0,90) 

title  (['Detail:  ',f]); 

%  asks  for  polarization  number  through  function  "filt" 
filt; 

but 3 =f indob j (gcbf , 'Tag' , 'profile')  ; 
set (but3, 'Enable' , 'on' )  ; 

case  'profile'  %  calculates  the  profile 

[x,y]=size(detail); 

prof ile=detail ( round (x/2 ) , :)  %%%integers 
%  transforms  digital  level  in  isothermal  units 
%  profile=prof ile*range/4095+level-range/2 ; 

%  transforms  digital  level  in  isothermal  units 
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for  z=l : length (profile) 

profile (z) =dl2temp (prof ile (z) , level, range, p) ; 
end 


plot (profile) 
grid 


%  enables  button  " 

peak " 

butl=f indob j (gcbf, 

'  Tag '  , 

' peak ' ) ; 

set (butl , ' Enable ' , 

'  on '  )  ; 

but3=f indob j (gcbf. 

'  Tag '  , 

'profile' )  ; 

set (but3 , ' Enable ' , 

'off') ; 

case  'peak'  %allows  selection  of  peaks  and  valleys 

title (' select  the  four  apparent  peaks,  using  the  mouse') 
[px,py] =ginput (4) ; 

title (' select  the  three  apparent  valleys,  using  the 
mouse' ) 

[vx, vy] =ginput (3 ) ; 
peak=mean (py) ; 
valley=mean (vy) ; 

text (mean (px) , peak, [ 'mean  peak  is: 

' , num2str (peak) , 1 °C ' ] ) ; 

text (mean (px) , valley, [ 'mean  valley  is: 

’ , num2str (valley) , ' °C ' ] ) ; 

title ([ 'Horizontal  profile  of  the  output  -  file:  ’,f]) 
butl=findobj (gcbf, 'Tag' , ’peak' ) ; 
set (butl, ' Enable ' , ' of f ’ ) ; 

%  enables  button  "fourier" 
but2=findobj (gcbf, 'Tag' , ' fourier' ) ; 
set (but2 , ' Enable ' , ' on ' ) ; 


case  'three'  %  plots  detail  in  3D 
surfl (detail) 
shading  interp; 
colormap (pink) 

title(['3D  view  of  the  output  -  file:  ’,f]) 

case  'fourier'  %  allows  frequency  domain  analysis 
minv=min(vy) ; 

shadow=f ind(profile/minv>=l) ; 
mask=zeros (1, length (prof ile) ) ; 
mask ( shadow) =1 ; 
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prof ilel=prof ile . *mask; 
mask (shadow) = (peak+ valley) / 2 ; 
prof ile2=prof i lei -mask; 


NET=0 . 12 ; 

y=abs (loglO (length (prof ile2) ) /loglO (2) )  ; 

F=f f t (prof ile2 ) / (2Ay) ; 

F=f ftshif t (F) ; 

mu= (-length (prof ile2) /2 : length (prof ile2 ) 12- 
1)  •“l/length(profile2)  ; 

plot (mu, abs (F) ) 
hold  on 

plot (mu,min (abs (F) ) *6*ones (length (prof ile2) ) ,'r') 
hold  off 

axis ( [min (mu)  max(mu)  0  max (max ( abs (F ) ) , min ( abs ( F ) ) * 6 ) ] ) 
but2=f indobj (gcbf, 'Tag' , 'fourier' )  ; 
set (but2 , ' Enable ' , ' off ' ) ; 

end 


% 

% 

% 

% 

% 

% 

% 


********* 


t************************************************* 


********* 

********* 

Type 

Name 

Function 


Naval  Postgraduate  School  -  CA 
Thesis  Research 
Function 
readgui 

Graphical  User  Interface  (GUI) 
function  readptw2.m  that  reads 
generated  by  PTRwin 


equivalent  to 
* .ptw  files 


%  Date  :  Jul  /  1999 

%  Version  :  1.0 

%  Author  :  Edson  Fernando  da  Costa  Guimaraes ,  Capt.  (BAF) 

%  Routine  Information : 

%  1)  loads  figure  from  file 

%  2)  gets  a  detail 

%  3)  calculates  the  profile 

%  4)  allows  selection  of  peaks  and  valleys 

%  5)  allows  frequency  domain  analysis 

%  6)  plots  detail  in  3D 

%********************************************************** 
$^*  **************************************  *  ■*•*■*  *  *  ************  * 


function  fig  =  readgui ( ) 
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%  This  is  the  machine-generated  representation  of  a  Handle 
%  Graphics  object  and  its  children.  Note  that  handle  values 
%  may  change  when  these  objects  are  re-created.  This  may 
%  cause  problems  with  any  callbacks  written  to  depend  on 
%  the  value  of  the  handle  at  the  time  the  object  was  saved. 
%  To  reopen  this  object,  just  type  the  name  of  the  M-file 
%  at  the  MATLAB  prompt.  The  M-file  and  its  associated  MAT 
%  -file  must  be  on  your  path. 

load  readgui 

hO  =  figure ( 'Units 'points ' ,  ... 

'Color' ,[0.80.80.8],  ... 

'Colormap' ,mat0,  ... 

' PointerShapeCData' ,matl,  . . . 

' Position' ,mat2 ,  ... 

' Render er' , 'zbuffer' ,  ... 

' RendererMode ' , 'manual' ,  ... 

'  Tag '  ,  '  Figl.' )  ; 

hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units ', 'points ' ,  ... 

' BackgroundColor ' , [0.752941176470588  0.752941176470588 
0.752941176470588],  ... 

' FontSize ' , 28 ,  ... 

'  Foregrounded  or '  ,  [0 . 3  0.30.3],... 

' ListboxTop' , 0 ,  ... 

'Position', [60.82758620689656  261.9310344827587 
325.8620689655173  32.89655172413794],  ... 

' String' , 'Welcome  to  Readptw.m' ,  ... 

'Style' , 'text' ,  ... 

' Tag ' , ' StaticTextl ' ) ; 
hi  =  uicontrol (' Parent ', hO ,  ... 

'Units' , 'points' ,  . . . 

'Callback' , 'guil ( ' ' load' ' ) ' ,  ... 

'ListboxTop' ,0,  ... 

'Position' ,  [290.4827586206897  205.448275862069 
72.62068965517243  29.17241379310346],  ... 

'String' , 'Load  PTW' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 
hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units' , 'points'  ,  ... 

'Callback' , 'guil ( ' ' three' ' ) ' ,  ... 

' ListboxTop ', 0 ,  ... 

'Position',  [290.4827586206897  64.55172413793105 
72.62068965517243  29.17241379310346],  ... 
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'String' , '3D  View' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 
hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units ', 'points ' ,  ... 

'Callback' , 'guil ( ' 'profile' ' ) ' ,  ... 

'Enable' , 'off' ,  ... 

' ListboxTop ' , 0 ,  ... 

' Position' ,mat3 ,  ... 

'String' , 'Profile' ,  ... 

'Tag' , 'profile' )  ; 

hi  =  uicontrol ( 'Parent ' ,h0,  ... 

'Units' , 'points ' ,  ... 

'Callback' , 'fclose( ' 'all' ' ) ; close (gcbf) ' ,  ... 

' ListboxTop' , 0 ,  ... 

'Position',  [290.4827586206897  19.24137931034483 
72.62068965517243  29.17241379310346],  ... 

' String' , 'Exit' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 

hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units' , 'points' ,  ... 

'Callback' , 'guil ( ' 'detail'  ' )  ' ,  ... 

' ListboxTop ' , 0 ,  ... 

'Position', [290.4827586206897  160.1379310344828 
72.62068965517243  29.17241379310346],  ... 

' String' , 'Detail' ,  ... 

' Tag ' , ' Pushbuttonl ' ) ; 
hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units' , 'points' ,  ... 

'Callback' , 'guil  (  "peak"  )' ,  ... 

'Enable' , 'off'  ,  ... 

' ListboxTop ', 0 ,  ... 

'Position', [39.10344827586208  16.44827586206897 
93.10344827586209  16.75862068965517],  ... 

' String' ,' Peak  &  Valley',  ... 

' Tag ' , ' peak ' ) ; 

hi  =  uicontrol ( 'Parent' ,h0,  ... 

'Units' , 'points' ,  ... 

'Callback' , 'guil( ' 'fourier' ' )  ' ,  ... 

'Enable' , 'off '  ,  ... 

' ListboxTop ' , 0 ,  ... 

'Position',  [157.0344827586207  16.44827586206897 
93.10344827586209  16.75862068965517],  ... 

'String' , 'Fourier  Transform' ,  ... 

' Tag ' , ' fourier ' ) ; 
hi  =  axes ( ' Parent ' , hO ,  ... 
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'Units' , 'pixels' ,  ... 

' Box ' , ' on ' ,  ... 

' CameraUpVector '  ,  [0  10],  ... 

'CameraUpVectorMode' , 'manual ' ,  ... 

' Color' ,[111],  ... 

'ColorOrder ' ,mat4,  ... 

' Position' , [66  106  337  271],  ... 

' Tag' , ' Axesl ' ,  ... 

' XColor ' ,[000],  ... 

' YColor ' ,[000],  ... 

' ZColor ' , [0  0  0] ) ; 
h2  =  line (' Parent ', hi,  ... 

'Color' ,[001],  ... 

'Tag'  ,  '  Axesl-Linel '  ,  ... 

' XData ' , 0 ,  ... 

' YData ' , 0 ) ; 

h2  =  text ( ' Parent ' , hi ,  ... 

'Color' ,[000],  ... 

'HandleVisibility ' , 'off ' ,  ... 

' HorizontalAlignment ' , 'center' ,  ... 

'Position',  [-0.005952380952380931  -1.207407407407407 

17.32050807568877],  ... 

'Tag' , 'Axes lText4 ' ,  ... 

' VerticalAlignment ' , ' cap' ) ; 
set  ( get  (h.2  ,  '  Parent '  )  ,  '  XLabel '  ,  h2 )  ; 
h2  =  text ( ' Parent ' , hi ,  ... 

'Color' , [0  0  0] ,  ... 

' HandleVisibility ',' of f ' ,  ... 

'HorizontalAlignment' , 'center' ,  ... 

'Position', [-1.238095238095238  -0.007407407407407085 

17.32050807568877],  ... 

'Rotation' ,90,  ... 

' Tag' , 'Axes lText3 ' ,  ... 

'VerticalAlignment' , 'baseline') ; 
set (get (h2, 'Parent' ) , 'YLabel' ,h2) ; 
h2  =  text ( ' Parent ' , hi ,  ... 

' Color' , [0  0  0] ,  ... 

'HandleVisibility' , 'off ' ,  ... 

'HorizontalAlignment' , 'right' ,  ... 

'Position',  [-1.392857142857143  1.940740740740741 

17.32050807568877],  ... 

'Tag' , 'Axes lText2 ' ,  ... 

'Visible ' , ' of f ' ) ; 

set (get (h2 , ' Parent ' ) , ' ZLabel '  ,  h2) ; 
h.2  =  text  ( '  Parent ' ,  hi ,  ... 
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'Color' ,[000],  ... 

' HandleVisibility ' , 'off' ,  ... 

'HorizontalAlignment ' , 'center'  , 
' Position' ,mat5 ,  ... 

'Tag' , 'AxeslTextl' ,  ... 

'VerticalAlignment ' , 'bottom' )  ; 
set (get (h2, 'Parent' ) , 'Title' ,h2) ; 
if  nargout  >  0,  fig  =  h0;  end 
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APPENDIX  B.  LABORATORY  SETUP 


In  order  to  measure  the  Minimum  Resolvable  Temperature 
Difference,  the  electro-optics  laboratory  was  prepared  to 
be  used  with  the  setup  shown  in  Figure  B.l  and  described  in 
this  Appendix. 


AGA  780 


1 


2 


3 


Figure  B.l  -  Laboratory  Setup  for  Measurement  of  MRTD. 

1.  Autotransformer  Variac®  type  WI0MT3  (General  Radio 
Company) ,  with  output  voltage  variable  from  0  to 
140V  (10  Amps) .  This  equipment  was  used  to  heat  the 
background  plate  by  varying  the  voltage  applied  to 
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a  resistance  installed  on  the  backside  of  the  plate 
(4)  . 

2.  Digital  Thermometer  model  2190A  (Fluke) .  This 
equipment  displays  either  Centigrade  or  Fahrenheit 
values,  both  with  a  precision  of  1  unit.  All  the 
measurements  were  taken  in  degrees  Fahrenheit  in 
order  to  improve  the  precision  (1°F~0 . 555°C) .  This 
thermometer  was  used  to  display  the  temperature  of 
the  back  plate  (4)  corresponding  to  the  background. 
The  equipment  precision  is  not  satisfactory  for 
MRTD  measurements.  Equipment  precision  should  be 
improved  if  further  measurements  are  to  be  taken. 

3.  Same  as  <2>,  but  used  to  display  the  temperature  of 
the  front  plate  (5) . 

4 .  Back  plate .  This  plate  represents  the  background 
with  a  emissivity  of  0.9. 

5.  Front  plate.  This  plate  represents  the  target  (kept 
at  room  temperature)  and  contains  the  standard  MRTD 
pattern  (7:1  aspect  ratio)  as  shown  in  Figure  B.2. 

6.  AGA  Thermovision  780.  It  is  a  two  channel  infrared 
scanning  system  with  single  detector.  It  operates 
in  both  MWIR  (3-5|im)  and  LWIR  (8-12(Jm)  [Ref.  4]  .  For 
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this  thesis,  only  the  LWIR  band,  with  the  HgCdTe 
detector,  was  used. 


7*a 


Figure  B.2  -  MRTD  Pattern. 

7.  Black  and  white  monitor  chassis.  It  was  used  as  an 
interface  between  the  camera  (6)  and  the  computer. 
Although  it  is  the  original  display  of  the  AGA  780, 
most  of  its  functions  were  performed  by  the  NEC 
computer  (8)  ,  which  provides  an  extended  set  of 
resources  for  recording  and  analyzing  the  images . 

8.  NEC  laptop  computer  VERSA  P75  (Pentium®  75  MHz) 
installed  with  PTRWIN®  v3.16  software.  This 
software  was  used  for  acquisition  and  analysis  of 
the  images  received  from  the  AGA  780.  Using  a  12- 
bit  A/D  converter,  it  converts  detector  signal 
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intensity  into  4096  discrete  digital  levels. 
12] 


[Ref . 
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APPENDIX  C.  SUMMARY  OF  THE  IMAGES 


Tables  C.l  to  C.3  show  a  summary  of  the  most  important 
data  extracted  from  the  images  collected  in  this 
experiment.  The  fields  in  this  table  are:  file  name  (File), 
spatial  frequency  (Freq) ,  background  temperature  (Tb 
converted  to  °C) ,  target  temperature  (Tt) ,  thermal  level 
(TL)  ,  thermal'  range  (TR)  ,  temperature  of  the  peak  (Tp)  , 
temperature  of  the  valley  (Tv) ,  MRTD,  calibrated  MRTD 
(MRTDcal) ,  polarizer  (p) ,  and  date  of  collection  (Date) . 


File 

Freq 

Tb 

Tt 

TL 

TR 

Tp 

Tv 

MRTD 

MRDTcal 

P 

Date 

A__lun00  .ptw 

0.43 

20 

52 

2 

0.54 

0 

10 /may 

A__lun01  .ptw 

0.45 

21.7 

52 

2 

2.37 

0 

10 /may 

A_lun02 .ptw 

0.47 

20.2 

22.3 

52 

2 

2.10 

4.28 

0 

10 /may 

A_JLun03 .ptw 

0.5 

20.2 

25.3 

54 

2 

25.79 

25.69 

5.10 

5.01 

0 

10 /may 

A_lun04 .ptw 

0.52 

21 

32.6 

56 

5 

29.41 

29.25 

11.60 

7.06 

0 

10 /may 

A_f 00 .ptw 

0.1 

24.3 

24.2 

51 

2 

20.15 

20.35 

0.05 

0 

2 5 /may 

A_f 02 .ptw 

0.2 

24.3 

24.6 

85 

2 

0.17 

0 

2 5 /may 

A_f 03 .ptw 

0.3 

24.3 

24.94 

51 

0.64 

0 

2 5 /may 

A_f 04 .ptw 

0.4 

2 

20.97 

0 

2  5 /may 

A_f 05 .ptw 

0.5 

24.44 

31.61 

52 

4.86 

0 

A— f 07 .ptw 

0.52 

24.61 

38.33 

56 

5 

10.13 

0 

A_unp00 .ptw 

0.3 

22.3 

54 

2 

26.26 

25.89 

1.90 

0.51 

0 

06 /may 

A_unp01 .ptw 

0.4 

22.2 

24.2 

54 

2 

0 

06 /may 

A_unp04 .ptw 

0.48 

23.5 

41 

59 

20 

35.00 

0 

06 /may 

A_unp06 .ptw 

0.43 

23.2 

5 

32.46 

32.13 

10.50 

3.20 

0 

06 /may 

A2__un01  .ptw 

0.1 

23.39 

23.33 

55 

0.05 

0 

11 /may 

A2_un02 .ptw 

0.2 

23.33 

23.61 

54 

2 

24.86 

24.70 

0.28 

0.17 

u 

A2_un03 .ptw 

0.3 

23.33 

24 

54 

24.95 

0.67 

0.62 

0 

11 /may 

A2_un04 .ptw 

0.4 

25.95 

1.28 

1.17 

0 

11 /may 

A2_un05 .ptw 

0.5 

23.72 

9.07 

0 

11 /may 

Table  C.l  -  Summary  of  Data  for  the  Unpolarized  Case. 
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File 

Freq 

Tb 

Tt 

TL 

TR 

TP 

Tv 

A_fp00 .ptw 

0.1 

23.61 

24.17 

49 

2 

21.35 

20.77 

A_fp01 .ptw 

0.2 

23.61 

24.33 

49 

2 

21.65 

21.25 

A_fp02 .ptw 

0.3 

23.61 

24.78 

50 

2 

25.63 

25.48 

A_fp03 .ptw 

0.4 

23.72 

25.78 

51 

25.94 

A_fp04 .ptw 

0.43 

23.72 

27.39 

2 

28.17 

27.94 

A_fp05 .ptw 

0.46 

23.78 

30.22 

52 

2 

30.13 

29.89 

A_fp06 .ptw 

0.5 

23.89 

37 . 61 

53 

2 

33.31 

33.07 

A  fp07.ptw 

0.52 

24.61 

65.56 

58 

5 

49.60 

49.42 

MRTD  I  MRDTcal 


0.10 


Table  C 


0.46|  23.78  30.22  52  2  30.13  29.89  6.44  2.67  4 

0.5  23.89  37.61  53  2  33.31~33.07  13.72  5.74  4 

0.52  24.61  65.56  58  5  49.60  49.42  40.95  23.06 _ 4 

,2  —  Summary  of  Data  (Vertically  Polarized  Case) . 


Date 


2 6 /may 


2  6 /may 


2  6 /may 


2  6 /may 


2  6 /may 


2  6 /may 
2 6 /may 
2  6 /may 


File 

Freq 

Tb 

Tt 

TL 

TR 

Tp 

Tv 

MRTD 

MRDTcal 

P 

Date 

A5_fp00 .ptw 

0.1 

16.94 

17.05 

9.06 

8.63 

0.11 

0.03 

5 

08/ jun 

0.2 

16.94 

17.28 

46 

2 

9.34 

9.12 

0.34 

0.16 

5 

08/ jun 

A5_fp02 .ptw 

16.94" 

47 

2 

13.12 

12.94 

0.84 

0.45 

m 

08/ jun 

A5__fp03  .ptw 

16.94 

19.39 

47 

2 

14.51 

14.37 

2.45 

1.76 

m 

08/ jun 

A5_fp04  .ptw 

0.43 

16.94 

2 

17.57 

17.34 

3.56 

1.53 

m 

A5_fp05 .ptw 

0.46 

16.94 

23.5 

49 

2 

21.08 

20.89 

6.56 

3.43 

5 

08 /jun 

A5_fp06 .ptw 

0.5 

17.11 

31.28 

49 

2 

22.92 

22.67 

14.17 

5.79 

5 

08/ jun 

A5_fp07 .ptw 

0.52 

17.11 

39.44 

51 

2 

27.27 

27.03 

22.33 

9.45 

5 

08/ jun 

Table  C.3  -  Summary  of  Data  (Horizontally  Polarized  Case) . 


Table  C.4  shows  the  summary  of  the  images  taken  for 
MTF  (Modulation  Transfer  Function)  estimation.  The  curve 
fitting  process  was  carried  out  using  the  program  mtffit.m 
(Appendix  A) . 


File 

Peak  (DL) 

Valley  (DL) 

Filter 

Mtf  00 .ptw 

4083.9416 

0.0000 

0 

MtfOl .ptw 

4083.9416 

0.0000 

0 

Mtf 02 .ptw 

3306.5693 

9.7324 

0 

Mtf 03. ptw 

2380.1703 

431.4680 

0 

Mtf 04 .ptw 

2753.6496 

1346.7153 

0 

Mtf 05 .ptw 

4083.9416 

4091.2409 

0 

Mtf 06 .ptw 

4083.9416 

4083.9416 

0 

Table  C.4  -  MTF  Data  Summary. 
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APPENDIX  D.  MEAN  RADIANCE  (8-12|im) 


T 

<N> 

nr 

<N> 

240 

2.861 12437e-004 

mm 

6.78699088e-004 

241 

2.9333 1304e-004 

iu 

6.91431324e-004 

242 

3.0067326  le-004 

282 

7.043 14046e-004 

243 

3.081 39 196e-004 

283 

7.17347725e-004 

244 

3.15729983e-004 

284 

7.30532826e-004 

245 

3.23446490e-004 

285 

7.43869799e-004 

246 

3.31289573e-004 

286 

7.57359088e-004 

247 

3.39260078e-004 

287 

7.71001 125e-004 

248 

3.47358840e-004 

288 

7.84796332e-004 

249 

3.55586687e-004 

7.98745 122e-004 

250 

3.63944432e-004 

8.12847900e-004 

251 

3.72432883e-004 

|  291 

8.271 0505 8e-004 

252 

3.81052833e-004 

wwm 

8.41516980e-004 

253 

3.89805068e-004 

mm 

8.5608404  le-004 

254 

3.9869036  le-004 

294 

8.70806605e-004 

255 

4.07709477e-004 

295 

8.85685028e-004 

4.16863169e-004 

9.00719655e-004 

■ 

4.26152180e-004 

wsm 

9.15910825e-004 

258 

298 

9.31258863e-004 

259 

4.45139079e-004 

9.46764089e-004 

260 

4.54838400e-004 

mm 

9.6242681  le-004 

261 

4.64675906e-004 

301 

9.78247330e-004 

262 

4.74652289e-004 

302 

9.94225936e-004 

263 

4.84768227e-004 

303 

1.01036291e-003 

264 

4.9502439  le-004 

304 

1. 02665 853e-003 

265 

5.05421439e-004 

305 

1.0431 1306e-003 

266 

5.15960020e-004 

306 

1.059726756-003 

267 

5.26640770e-004 

307 

1.07649985e-003 

268 

5.374643 17e-004 

308 

1.09343260e-003 

269 

5.48431278e-004 

309 

1.11052522e-003 

270 

5.59542259e-004 

310 

1.12777795e-003 

271 

5.70797856e-004  I 

311 

1.14519099e-003 

272 

5.82198655e-004 

312 

1.16276456e-003 

273 

313 

1.18049883e-003 

274 

6.05438 143e-004 

ItlW 

1.19839401e-003 

275 

6.17277952e-004 

Hmi 

1.21645027e-003 

276 

6.29265200e-004 

316  | 

1.23466779e-003 

277 

6.41400419e-004 

BU 

1.25304673e-003 

278 

6.53684133e-004 

mm 

1.27158725e-003 

279 

6.661 16855e-004 

319 

1.29028950e-003 

320 

1.30915361e-003 

Table  D.l  -  Mean  Radiance. 
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